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EXAMINATION  AND  EXPERIMENTAL  EVALUATION 
OF  A  MATHEMATICAL  ANALYSIS  OF  PULSE  TUBE  REFRIGERATION  ' 

by 
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Eugene  Edward  Fitzpatrick,  Lieutenant,  United  States  Navy 
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degree  of  Master  of  Science  in  Mechanical  Engineering  and  the  profes- 
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ABSTRACT 

An  examination  of  a  mathematical  analysis  of  Pulse  Tube  refrig- 
eration as  proposed  by  Rea  was  carried  out.  Rea  originally  investi- 
gated the  performance  of  thermal  regenerators  subjected  to  rapid 
pressure  and  flow  cycling  and  found  that  this  class  of  regenerators 
is  characterized  by  two  basic  parameters; 

(1)  The  ratio  of  the  convective  heat  transfer  within  the  regenerator 
to  the  heat  generated  by  compression, 

(2)  The  ratio  of  the  rate  of  mass  storage  within  the  void  volume  of 
the  regenerator  to  the  through-flow. 

Rea's  proposal  was  to  treat  the  Pulse  Tube  as  a  very  efficient  regen- 
erator in  series  with  a  very  inefficient  regenerator,  and,  thus, 
characterize  the  performance  of  the  Pulse  Tube  by  the  same  two  para- 
meters. 

The  examination  of  Rea's  proposal  indicates  that,  in  addition 
to  the  two  parameters  which  describe  the  performance  of  an  efficient 
regenerator,  there  is,  for  the  case  of  the  inefficient  regenerator,  a 
third  parameter  which  depends  upon  the  average  pressure  and  the  rate 
at  which  the  pressure  and  mass  flow  cycle. 

An  experimental  evaluation  of  Rea's  proposal  modified  to  include 
the  third  parameter  was  performed  on  a  Pulse  Tube  consisting  of  a 
spherical  matrix  regenerator  in  series  with  a  void  tube.  There  was 
good  agreement  observed  between  theoretical  and  experimental  results* 

Thesis  Supervisor:  Joseph  L.  Smith,  Jr. 

Title:  Associate  Professor  of  Mechanical  Engineering 
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NOMENCLATURE 
A     -    cross-sect3.onal  area  of  the  regenerator  or  pulse  tube, 

«2. 

A0    -    cross-sectional  area  of  the  regenerator  or  pulse  tube  open 
to  flow  (A0  =  £  A),  ft2. 

Arj,  -    heat  transfer  surface  area  per  foot  of  length  of  regenerator 

or  pulse  tube,  ft  /ft. 

B     -     dimensionless  variable  defined  by  equation  (A.U8). 

C     -     dimensionless  constant  defined  by  equation  (2.31). 

C0    -    integration  constant  (C0  «  01To   -  C1/Te). 

Cn     -     constant  defined  by  equation  (2.19),  lbm  -  °F/hr. 

cm    -     specific  heat  of  regenerator  matrix  or  pulse  tube  wall, 
BTU/lb™  -  °F. 

c     -     constant  pressure  specific  heat  of  the  gas,  BTU/lbm  -  F. 

c     -     constant  volume  specific  heat  of  the  gas,  BTU/lbm  -  °F. 

h     -    gas  enthalpy  at  a  position  x  along  the  regenerator  or  pulse 
tube,  BTU/lbm. 

h,hm   -    heat  transfer  coefficient  at  position  x  along  the  regenera- 
tor or  pulse  tube,  BTU/hr-ft2  -  °F. 

K     -     heat  transfer  coefficient,  BTU/hr-ft2-°F. 

k     -     thermal  conductivity,  BTU/hr-ft-°F0 

L     -     regenerator  or  pulse  tube  length,  ft. 

m     -    gas  mass  flow  at  position  x  along  regenerator  or  pulse 
tube,  lb  /hr. 

m     -    dimensionless  mass  flow  (m  =  m/m0). 
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mQ    -    mass  flow  at  x  z   0  in  regenerator  or  pulse  tube,  lbm/hr# 
M     -    matrix  mass  per  foot  of  regenerator  length  or  mass  per 

foot  of  pulse  tube  length,  lbm/ft. 
n     -    positive  exponent  whose  value  depends  on  regenerator 

matrix  (n  =  0  for  laminar  flow  pulse  tube). 

p     -     gas  pressure  in  regenerator  or  pulse  tube,  lbf/ft  • 

Q     -     heat  flux,  BTU/hr. 

R     -     gas  constant,  ft-lb„/lb  -  F. 

i   m 

t     -     time,  hours, 

T     -     temperature  of  the  regenerator  matrix  or  pulse  tube  wall 

at  position  x,  F„ 
T     -    temperature  of  the  gas  at  position  x  along  the  regenera- 

tor  or  pulse  tube,  (T. 
u     -    internal  energy  of  the  gas,  BTU/Lbm. 
V     -    volume,  ft3. 

V„    -     dead  volume  at  end  of  pulse  tube,  ft3. 
Vp^    -    volume  of  pulse  tube,  ft-^o 
Vb     -    void  volume  of  regenerator,  ft^» 

VT    -     total  volume  of  the  regenerator  or  pulse  tube,  ft3, 
w     -     dimensionless  distance  along  regenerator  or  pulse  tube 

defined  by  equation  (2,26). 
x     -     distance  along  regenerator  or  pulse  tube,  ft. 


Greek  Letters 

oc    -     dimensionless  Pulse  Tube  parameter  defined  by  equation  (3.18), 
8  -     heat  transfer  area  per  unit  total  volume  of  the  regenerator 

or  pulse  tube  (  f  =  ~)y   ft2/ft3, 

-x- 


£ 

$ 

9 

e 
r 
V 


dimensionless  regenerator  or  pulse  tube  void  fraction, 
dimensionless  ratio  of  specific  heats  (If  :  c  /c  ). 
dimensionless  pulse  tube  parameter  defined  by  equation 

(2.3*). 

gas  density,  lbjj/ft  • 

dimensionless  matrix  or  pulse  tube  wall  temperature 

(0=  T/T0). 

duration  of  cycle,  hours, 

duration  of  a  compression  period,  hours. 

duration  of  an  expansion  period,  hours, 

percentage  of  the  total  cycle  time  devoted  to  compression 

or  expansion. 


Subscripts 

avg    -  average o 

c     -  compression, 

C     -  conditions  at  cold  end  of  Pulse  Tube, 

cond   -  conduction, 

e     -  expansion, 

eff    -  effective, 

exper  -  experimental, 

g     -  gas. 

heater  -  used  with  Q  to  indicate  heater  input  to  Pulse  Tube. 

H     -  conditions  at  hot  end  of  pulse  tube. 

losses  -  used  with  Q  to  indicate  conduction  and  radiation  losses 

in  Pulse  Tube, 

m     -  refers  to  regenerator  matrix  or  pulse  tube  wall. 
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max  -  refers  to  maximum  value  of  the  variable  w  or  p. 

min  -  refers  to  minimum  value  of  the  variable  p. 

o  -  pulse  tube  or  regenerator  conditions  at  x  -   w  s  0, 

pt  -  pulse  tube0 

rad  -  radiation. 

o 

R  -  regenerator;  also  refers  to  refrigeration  (Q^)« 

theor  -  theoretical, 

¥  -  conditions  at  warm  end  of  regenerator. 
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CHAPTER  I 


bitroduct: 


1 

A  recent  paper  by  Gifford  and  Longsworth  (l)"*  describes  an  inter- 
esting thermodynamic  effect  which  frequently  is  ignored  in  practice, 
but  '.  hich  has  a  potentially  useful  applicationo  As  Gifford  and 

FIGURE  I 

SCHEMATIC  OF  THE  PULSE  TUBE 
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Longsworth  indicate  (2),  it  is  not  generally  recognized  that  the  cycling 
of  a  gas  in  a  closed  volume  can  establish  large  temperature  gradients 
within  such  a  volume.  Gifford  and  Longsworth,  by  suitably  arranging  a 
thermal  regenerator,  heat  exchangers,  and  a  void  volume,  have  created  a 


Numbers  in  (  )  refer  to  References, 
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i(  3  '  hich  preserves  these  temperature  gradients  in  an  essentially 
tic  state,  even  though  the  gas  is  rapid  cyclir  \,     They 

/e  called  this  device  a  Pulse  Tube.  It  is  illustrated  sch   iically 
i  .  Fig   ■  I. 

^ruction  and  operation  of  the  device  are  simple.  This 
simplicity  when  coupled  with  the  promising  experimental  results  obtained 

is  far,  indicates  that  the  Pulse  Tube  should  find  practical  applica- 
tion in  the  near  future,  perhaps  as  a  small  cooler  in  electronic  equip- 

...  Gifford  and  Longsworth  (2)  report  refrigeration  across  temperature 
intervals  as  large  as  209°  F.  with  a  single-stage  system.  The  units  may 
also  be  operated  in  multi-stages,  so  larger  temperature  differences  are 
entirely  feasible. 

3  Pulse  Tube  consists  of  a  regenerator  connected  in  series  with 
a  void  tube,  which  is  the  pulse  tube.*""'"  The  void  volume  of  the  regen- 
erator is  Vp#  The  volume  of  the  pulse  tube  is  V  . .  Gas  enters  and 
rv  c  pt, 

exits  the  device  through  the  warm  end  of  the  regenerator,  which  is  at 

iiperature  T...  The  volume  Vv,  which  is  connected  to  the  end  of  the 
pulse  tube,  is  maintained  at  a  constant  temperature,  T,..  This  may  be 
accomplished  by  circulating  cooling  water  around  the  volume  Vw. 

The  device  operates  in  the  following  manner:  the  inlet  valve 
opens  admitting  high  pressure  gas  into  the  regenerator,  thereby 
inducing  a  mass  flow  into  the  regenerator  and  pulse  tube. 


Pulse  Tube  capitalized  r.     to  the  entire  device  (regenerator 
plus  pulse  tube);  pulse  tube  without  capital  letters  refers  to 

.  void  tube.  Wherever  the  void  tube  is  intended,  "cut  gram- 
matically capital  letters  are  required,  the  following  notation 
will  be  used:   "pulse  tube". 
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CU3  r  vtly,  there  is  a  continuous  pressure  rise  throughout  the  Pulse 

some  predetermined  maximum  tires  sure,  p   ,  within  the  device, 

1,1  nay:3 

inlet  valve  closes  and  the  exit  valve  opens.  This  allows  the  gas 
hin  the  Pulse  Tube  to  flow  back  out  through  the  warm  end  of  the  regen- 
erator, 

an  the  inlet  valve  is  open,  the  gas  in  the  pulse  tube  is  heated 
as  it  is  compressed  by  the  entering  gas.  Some  of  this  heated  gas  flows 
into  the  volume  VT ,  where  it  is  cooled  by  the  cooling  medium  to  the  tem- 
perature T„.  When  the  exit  valve  opens,  the  gas  is  cooled  further  by 
expansion  as  it  flows  back  out  of  the  volume  V  ,  Some  of  this  cool  gas 
must  precool  the  entering  warm  gas  during  the  next  compression  half- 
cycle.  The  remainder  of  the  cool  gas  is  available  for  refrigeration, 

:rigeration  becomes  available,  once  a  steady  cyclic  state  is  reached, 
at  the  connecting  tubing  between  the  regenerator  and  the  pulse  tube* 
This  refrigeration  load  is  labeled  Q  in  Figure  I  and  is  available  at 
the  temperature  Tp. 

Previous  Work 

The  first  analysis  of  the  pulse  tube  was  given  by  Gifford  and 
Longsworth  in  their  first  and  second  papers  (1,2),  This  analysis  did 
not  predict  compatible  pulse  tube  capacity  results  with  those  they  had 
obtained  experimentally,  Gifford  and  Longsworth  recognized  that  their 
original  analysis  of  the  pulse  tube  was  inadequate,  Tneir  original 
analysis  neglected  any  heat  exchange  between  the  gas  and  the  pulse  tube 
wall.  They  discovered  that  this  heat  exchange  could  contribute  a  sub- 
stantial amount  to  the  refrigeration  effect.  In  other  words  the  pulse 
tube  itself  acts  as  an  inefficient  regenerator.  The  third  paper  by 
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Gifford  and  Longsworth  (3)  discussed  the  aspects  of  this  heat  exchange 
ween  the  gas  and  the  pulse  tube  vail  in  detail.  However,  no  complete 
lysis  including  this  effect  was  presented.  The  problem  was  primarily 
.  of  determining  the  appropriate  temperature  difference  to  use  in 

FIGURE  II 
SCHEMATIC  0?  REA'S  2-PART  MODEL  FOR  A  REGENERATOR  (U) 


Compressor 


Inlet  Valve 


Regenerator 
Void  Volume  r  7 


R 


order  to  calculate  correctly  the  heat  transferred  in  the  process  under 
consideration. 

Concurrently  with  the  work  of  Gifford  and  Longsworth,  Rea  was 
studying  the  performance  of  thermal  regenerators  subjected  to  rapid 
pressure  and  flow  cycling  (U).  Rea's  objective  was  to  develop  a  theory 
for  predicting  the  axial  temperature  distribution,  mass  flow  rate 
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iion,  and  effectiveness  of  thermal  regenerators  operated  under 
the  above  conditions.  In  developing  his  theory  Rea  established  his  so- 
called  2-Part  Model  (h),   which  is  patterned  after  the  Gifford-McMahon 
cycle  (£)•  A  schematic  of  Rea's  2-Part  Model  is  shown  in  Figure  II. 

By  making  energy  and  mass  balances  for  the  gas  and  the  matrix  in 
an  elemental  control  volume ,  Rea  was  able  to  write  the  governing  differ- 
ential equations  for  his  problem.  After  applying  appropriate  simplifi- 
cations, Rea  reduced  his  problem  to  a  system  of  two  differential  equa- 
tions, which  he  solved  numerically  with  the  aid  of  a  computer.  A  complete 
description  of  Rea' s  2-Part  Model  and  the  development  of  the  governing 
differential  equations  are  presented  in  Appendix  A. 

When  the  work  of  Giff ord  and  Longsworth  became  known  to  him,  Rea 
recognized  that  the  Pulse  Tube  operates  in  the  same  manner  as  that 
described  by  his  2-Part  Model.  Rea  concluded  that  since  a  pulse  tube 
was  nothing  more  than  an  inefficient  regenerator,  his  theory  should 
apply  equally  well  to  the  entire  Pulse  Tube  as  to  a  simple  regenerator. 
Thus,  Rea  extended  his  theory  to  include  the  pulse  tube. 

Rea  tested  his  theory  for  the  regenerator  and  got  fairly  accurate 
results.  Rea,  however,  did  not  test  his  theory  as  extended  to  include 
the  pulse  tube.  Thus,  it  is  the  purpose  of  this  thesis  to  determine 
whether  or  not  Rea's  Theory  correctly  describes  the  performance  of  the 
Pulse  Tube  by,  first,  examining  his  theory  and  analysis  for  the  entire 
Pulse  Tube  and  by,  second,  measuring  the  refrigeration  effect  of  an 
operating  Pulse  Tube. 
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CHAPTER  II 

EXAMINATION  OF  REA'S  TIIFORY  AS  APPLIED 

TO  A  REGEMFRA.TOR  AND  "pulse  tube" 

Rea's  2 -Part  Model  (U) 

When  faced  with  the  task  of  formulating  a  mathematical  analysis 
of  a  given  system,  the  most  important  step  is  the  selection  of  an 
appropriate  model o  Ideally  one  wants  a  model  that  is  amenable  to  a 
theoretical  treatment.  On  the  other  hand  one  does  not  want  his  model 
to  be  over-simplified,  for  then  the  problem  may  lose  its  significance. 
Rea  selected  as  his  model  for  the  regenerator,  and,  hence,  the  pulse 
tube,  a  model  that  satisfies  both  of  these  criteria.  This  "2-Part 
Model",  as  Rea  refers  to  it,  is  shown  schematically  in  Figure  II  of 
Chapter  I.  The  usefulness  of  this  model  is  proved  by  the  fairly 
accurate  results  which  are  obtained  by  its  use  and  by  the  fact  that 
certain  dimensionless  parameters  which  are  definite  aids  in  defining 
the  problem  evolve  from  its  use. 

The  operation  of  the  2-Part  Model  is  as  follows:  the  inlet 
valve  opens  allowing  high  pressure  gas  to  enter  at  the  top  of  the 
device.  During  this  half -cycle  the  pressure  within  the  device  con- 
tinually increases  with  the  mass  flow  rate  out  of  the  bottom  of  the 
device  equal  to  mQ.  After  a  predetermined  maximum  pressure,  p  max> 
is  reached  within  the  device,  the  inlet  valve  closes  and  the  exit 
valve  opens.  During  the  exhaust  half-cycle  the  process  described 
above  is  reversed,  with  mass  flow  rate  into  the  bottom  of  the 
device  equal  to  m  .  The  pressure  within  the  device  decreases  to 
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its  minimum  value,  Pj^,  during  this  half-cycle.  After  a  period  of 
time,  a  steady,  cyclic  condition  is  reached. 

In  examining  Rea's  theory  for  the  Pulse  Tube,  the  theory  for  the 
individual  components,  i.e.,  the  regenerator  and  pulse  tube,  will  be 
examined  separately. 

Pasic  Differential  Equations  for  the  Regenerator  and  "pulse  tube" 

The  governing  differential  equations  are  derived  in  detail  in 
Appendix  A  by  making  a  mass  balance  and  energy  balance  for  the  gas  and 
for  the  matrix  in  an  elemental  control  volume,  dx,  at  a  position  x 
within  the  regenerator  and  pulse  tube  respectively.  The  derivation  of 
the  regenerator  equations  is  entirely  attributable  to  Rea  (U).  The 
modification  to  Rea's  Theory  as  applied  to  a  pulse  tube  is  originalo 
The  basic  equations  for  the  regenerator  and  pulse  tube  are: 

Gas  Energy  Balance; 

cp  A  <;Tg)  ♦  hAT(Tg  -  t)  ♦  c^  a? 
0  r     at         * 


Matrix  Energy  Balancer 
Gas  Continuity: 


hA  (T  -T)=Mcm^  ,       (2#2) 


ax "    r  at\  y  "      r  \Tgat 


(2.3) 


-7- 


The  assumptions  necessarily  made  in  deriving  equations  (2,1), 
(2.2),  and  (2,3)  are: 

(1)  Zero  thermal  conductivity  in  the  direction  of  flow  for 
both  the  gas  and  the  matrix  and  infinite  thermal  con- 
ductivity in  the  matrix  transverse  to  the  flow  direction, 

(2)  The  pressure  drop  in  the  device  is  neglected, 

(3)  The  gas  obeys  the  perfect  gas  law, 

(k)     The  gas  transport  properties  and  the  matrix  properties 
are  constant, 

(3>)  The  cross  sectional  area  of  the  device  open  to  flow,  AQ, 
is  constant. 

In  addition  to  these  assumptions  it  is  necessary  to  assume  that 
the  heat  transfer  area  per  unit  volume  of  the  pulse  tube  is  fairly 
large,  in  order  that  the  pulse  tube  can  be  considered  as  a  regenera- 
tor. 

Boundary  conditions  for  equations  (2.1),  (2.2),  and  (2.3)  for 
the  regenerator  are: 

at  x  =  0:    m  =  m0(t) 

Tg  =  Tgo(t)  (2.1*) 

at  x  =  L:   Tg  =  TgW(t)  . 

Boundary  conditions  for  equations  (2.1),  (2.2),  and  (2.3)  for 
the  pulse  tube  are: 

at  x  =  0:    m  =  m0(t) 

Tg  =  Tg0(t)  (2.5) 


at  x  =  L:   Tg  =  TgC(t) 
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Note  that  these  boundary  conditions  are  identical  and  have  been 
presented  separately  only  to  show  that  the  position  x  -  L  for  the 
regenerator  is  the  warm  end  of  the  regenerator,  while  the  position 
x  =  L  for  the  pulse  tube  is  the  cold  end  of  the  pulse  tube,  when  the 
Pulse  Tube  is  operated  as  a  refrigerator. 

For  a  regenerator  (and  pulse  tube  in  this  case)  operating  in  the 
steady,  cyclic  manner  described  here,  certain  cyclic  conditions  must  be 
met.  At  a  given  position  x,  these  conditions  are  as  follows: 

<£hAT(Tg  -  T)dt  =  Mcm    *1  dt  =  0  (2.6) 

ot 

<^mdt  =  0  (2.7) 

$dt=0  (2.8) 

^(mTg)dt  =  C0   .  (2.9) 

Equation  (2.6)  says  that  over  a  complete  cycle  there  can  be  no 
net  heat  stored  by  the  matrix.  Equation  (2.7)  says  that  there  is  no 
net  mass  stored  over  a  complete  cycle.  Equation  (2.8)  states  that  the 
pressure  is  cyclic.  Equation  (2.9)  results  from  equations  (2.6)  and 
(2.8)  applied  to  a  cyclic  integral  of  equation  (2.1).  CQ  is  constant 
with  both  time  and  position.  The  cyclic  conditions,  equations  (2.6) 
to  (2.9),  are  general  in  nature  and  must  be  satisfied  by  any  steady, 
cyclic  solution  to  equations  (2.1),  (2.2),  and  (2.3). 

The  complete  cycle  will  now  be  broken  into  two  half -cycles  of 
equal  time  duration.  Also,  the  pressure  is  assumed  to  be  a  sawtooth 


function  of  time  such  that 


dp 
dt 


z  a  constant.  Other  approximations  and 
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assumptions  are  as  foil  oris: 

(1)  The  temperature  difference  between  the  gas  and  the  matrix 
at  any  point  within  the  regenerator  (or  pulse  tube)  is 
very  small.  This  enables  one  to  substitute  T  for  T  in  the 
final  continuity  relations. 

(2)  The  product  of  the  average  equals  the  average  of  the  product. 
This  approximation  is  useful  in  applying  the  cyclic  condi- 
tions, equations  (2,6)  through  (2.9) • 

(3)  For  the  regenerator  only, 


T  J   T  dt 


With  these  conditions  imposed,  the  final  set  of  equations  to  be 
solved  for  the  regenerator  are: 


-Cp  L.  (mTgc)  +  hAT(Tgc  -  T)  +  CyAp 


dx 


R 


dp 
dt 


Z  0       ,     (2.10) 


cp  L.  (mTge)  +  hAT(Tge  -  T)  -^o 
dx  T> 


dp 
dt 


=  0         ,     (2.11) 


• 

dm  _  A0 
dx  ""  RT 

dp 
dt 

t 

m(Tgc  ~  Tge)  =  cl 

T   +  T 
xgc  *  xge 

=  2 

T 

(2.12) 

(2.13) 
(2.11*) 


The  boundary  conditions  for  these  equations  are; 


at  x  z  0: 


at  x  :  L: 


T  z  T 


m  z  m 


T  = 


\r 


(2.15) 
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The  final  set  of  equations  to  be  solved  for  the  pulse  tube  are: 


■cp  J.  (mT     )  ♦  hAjT  ^  -  T)  +  Vo 
dx         gc 


Tx  gc 


R 


c_.A, 


cp  $g  (^ge)  +  hA^T^  -  T)  -  fv 


Tv   ge 


R 


dp 
dt 

dp 
dt 


—  -la 

dp 

dx       RT 

dt 

^V  "  Tge}  = 

Tgc  +  Tge 

=  2T 

1  -  2C    JL\  J2    i  I  Ha 
Uf-l/At     T    a 


=  c. 


:o     , 


:  0    , 


V 
KATL/ 


(2.16) 
(2.17) 

(2.18) 

(2.19) 
(2.20) 


The  boundary  conditions  for  these  equations  are: 


at  x  =  0: 


at  x  r  L: 


T  =  T 


m  -  itl 


T  =  Tr 


(2.21) 


Note  that  the  set  of  equations  for  the  regenerator  is  identical 
with  the  set  of  equations  for  the  pulse  tube  with  the  exception  that 
the  continuity  equations  (2.12)  and  (2.18)  are  different.  The  differ- 
ence between  these  two  equations  is  that,  in  the  case  of  the  regenera- 
tor, the  second  term  of  equation  (2.3)  can  be  shown  to  be  negligibly 
small  in  comparison  to  the  first  term,  whereas  in  the  case  of  the 
pulse  tube,  the  second  term  is  not  small  for  all  conditions  of 
operation.  Rea  assumed  that  the  continuity  equation  as  written  for 
the  regenerator  applied  equally  to  the  pulse  tube.  The  details  of 
the  relative  orders  of  magnitude  of  the  terms  of  equation  (2.3)  for 
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the  regenerator  and  pulse  tube  are  presented  in  Appendix  A, 

Equations  (2.10)  through  (2.11j)  are  five  equations  in  five  un- 
knovms  for  the  regenerator  and  can  be  reduced  to  two  equations  in  two 
unknowns  by  subtracting  (2.11)  from  (2,10),  and  then  by  substituting 
into  this  result  equations  (2.12),  (2.13),  and  (2.1U),  plus  the  com- 
monly known  relationship  R  =  cD  -  c^,.  The  result  of  this  manipulation 
is: 


i£  _  hAfl  +^o 


dx 


2cpm 


dp 
dt 


-  0 


(2.22) 


and, 


dm   A^ 
_  o 

dx   RT 


dp 
dt 


(2.23) 


The  three  boundary  conditions  (2.15)  are  required  since,  although  C^ 
in  equation  (2.22)  is  a  constant,  it  is  not  known  at  this  point. 

Idkexd.se,  equations  (2.16)  through  (2.20)  are  five  equations  in 
five  unknowns  for  the  pulse  tube  and  in  a  similar  fashion  to  that 
described  above  can  be  reduced  to  two  equations  in  two  unknowns. 
These  results  are: 


•  dT   hAmCT   A. 
m  —  -   L     ,  +  —2 

dx   2c  m    c 
P     P 


dp 
dt 


U-l/At  T  \m 


JT-1 


-  1 


=  0 


(2.2U) 
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and, 


dm  _  AQ 
dx   RT 


dp 
dt 


1  -  2C 


P  1  In, 


V 


ky-l  j  At  Tin  /  \  KATL 


(2.25) 


Again,  the  three  boundary  conditions  (2,21)  are  required,  since  C-j_  is 
unknown  at  this  point. 

Set-up  of  the  Differential  Equations  for  Computer  Solution 

The  following  dimensionless  variables  are  now  introduced: 

G  Z   T/T 


m  Z 


m/mo 


w 


=   V 


TqAqR 


dp 

dt 


(2.26) 


x 
L 


Also  it  is  assumed  that  the  heat  transfer  coefficient  can  be 
correlated  as: 


h  =  Knf1 
for  the  regenerator,  and  as: 

h  =  K 


(2.27) 


(2.28) 


for  the  pulse  tube,  which  is  assumed  to  be  operating  in  laminar  flow. 

When  the  new  variables  (2.26),  along  with  the  assumption  (2.27), 
are  introduced  into  the  regenerator  equations  (2.22)  and  (2.23),  the 
final  regenerator  equations  result  in: 


dO- 
dw 


m 


2-n 


(2.29) 


-13- 


and 


dm  m  1 
dw  ~  Q 


(2.30) 


where 


c  =  y-i  /katl  \  Cx       1 


V    /    2m, 


dt 


(2.31) 


The  boundary  conditions  for  the  regenerator  become: 


at  w  :  0:       m  =  1 

e=i. 

V 
at  w  =  w  _.      = 

dp 

dt 

V 

n0R 

:  0:1:6 


w 


(2.32) 


In  a  similar  fashion,  when  the  new  variables  (2.26),  along  with 
the  assumption  (2.28),  are  introduced  into  the  pulse  tube  equations 
(2.2U)  and  (2.25),  the  final  pulse  tube  equations  result  in: 


l  +  3) 

dQ-  £J J[ 

dw 


Jf-J. 

IT" 


(2.33) 


tor 


and 


~  =  i  I  1  -l£ 

^  0  I   0m 


(2.31;) 


where  C  is  defined  by  equation  (2.31)  and 


$« 


=  2 


JL.  R.L. 

KAyL  J  At  TQ 


(2.35) 
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The  boundary  conditions  for  the  pulse  tube  become: 
at  w  =  0:    m  =  1 

e=i 

0  0  o 

Method  of  Solution 

The  method  of  solution  for  the  regenerator  eauations  (2.29)  and 
(2.30)  is  to  treat  C  as  a  parameter.  This  makes  the  problem  an  initial 
value  problem,  since  both  ©  and  m  are  known  at  w  =  0,  To  obtain  a  set 
of  solutions,  one  must  know  the  ratio  of  the  specific  heats,  X    ,  and 
the  type  of  matrix  used  in  order  that  n  can  be  assigned  a  value.  The 
parameter  C  is  then  assigned  a  range  of  values,  and  the  equations  are 
numerically  integrated  from  w  =  0  to  w  =  wmax.  Solutions  for  n  =  ,£9 
and  tf  -   1,67  are  presented  in  graphical  form  in  Appendix  Ca 

The  method  of  solution  for  the  pulse  tube  equations  (2,33)  and 
(2 .3U)  is  similar  to  that  for  the  regenerator,  Now,  however,  there  is 
a  second  parameter,  Q>  ,  involved.  In  this  case  (|>  is  assigned  a  value, 
O   Z   C£>  ,  and  then  C  is  assigned  a  range  of  values  for  this  <J>  o  Again 
the  eouations  are  numerically  integrated  from  w  =  0  to  w  =  w^-^.  The 
result  is  a  family  of  curves  of  ©  versus  w  and  m  versus  w,  with  each 
curve  having  a  specific  value  of  C  and  each  graph  representing  one 
value  of  the  parameter  Q  ,  Solutions  for  various  values  of  <£  and 
#  =  I.67  are  presented  in  graphical  form  in  Appendix  C, 

To  use  the  curves  in  Appendix  C  one  must  know  wrnax  and  G^ 
for  the  case  of  the  regenerator.  In  the  case  of  the  pulse  tube, 
Wr   >  0  «,  and  (fe   must  be  known.  Then  entering  the  appropriate 
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graph  with  w  =  Wj^  and  G  =  @w,  or  ©  =  QQ,   a  value  of  C  can  be 
read.  The  axial  temperature  distribution  is  determined  by  following 
that  value  of  C  down  to  w  =  0.  Also,  once  C  is  known,  the  mass  flow 
into  the  top  of  the  regenerator  (or  the  bottom  of  the  pulse  tube,  as 
the  case  may  be)  is  determined  by  entering  the  graph  of  m  versus  w  at 
w  =  wmax  anc*  reading  off  the  value  of  m  corresponding  to  the  inter- 
section of  the  C  curve  and  the  w  =  w^^  line*  As  in  the  case  of  the 
temperature  distribution,  the  axial  mass  flow  distribution  is  deter- 
mined by  following  the  same  C  curve  from  w  =  wrcax  to  w  =  0.  This 
yields  the  mass  flow  distribution  from  the  warm  end  to  the  cold  end 
of  the  regenerator  or  from  the  cold  end  to  the  warm  end  of  the  pulse 
tube. 

From  this  point  on  equations  (2.29),  (2.30),  (2.33),  and  (2.3U), 
plus  the  boundary  conditions  equations  (2.32)  and  (2.36),  will  be 
referred  to  as  Rea's  Modified  Theory  for  the  Pulse  Tube. 
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CHAPTER  III 


ANALYSIS  OF  THE  PULSE  TUPE 


The  System 

The  Pulse  Tube  consists  of  a  regenerator  in  series  with  a  pulse 
tube,  as  shown  schematically  in  Figure  III,  This,  then,  is  the  system 
which  will  be  analyze d# 

FIGURE  III 
WORKING  SCHEMATIC  OF  THE  PULSE  TUBE 
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Analysis 

If  one  considers  the  control  volume  around  the  entire  Pulse  Tube 
as  shown  in  Figure  III  and  writes  an  energy  balance  for  one  cycle  of 
length  T',  the  result  is: 


Tccp  (Vgc)w  +  heater T  +  bosses r  "  ^ecp 


m  T 
e  ge 


W 


.  QHr=  o     ,  (3.1) 


where  the  subscript  W  refers  to  the  warm  end  of  the  regenerator.  In 
equation  (3»l)  the  first  terra  represents  the  enthalpy  carried  by  the 
gas  into  the  control  volume  at  the  warm  end  of  the  regenerator,-  the 
second  term  is  a  known  heat  quantity  introduced  into  the  control  volume, 
the  third  term  is  the  heat  associated  with  conduction  and  radiation  that 
is  absorbed  by  the  control  volume,  the  fourth  term  represents  the  enthalpy 
carried  by  the  gas  out  of  the  control  volume  at  the  warm  end  of  the  regen- 
erator, and  the  fifth  term  is  the  heat  removed  from  the  volume  V^# 

Now,  if  the  control  volume  around  the  pulse  tube  only  is  considered, 
one  can  write  an  energy  balance  for  one  cycle  of  length  T"  as  follows: 

T  c   mT  |   +Q1.     f-Tc   nl    -  Q„T  =  0  ,  (3.2) 
'  c  p  \   c  gc|c    losses     '  e  p  \   e  gejc    H      ' 

where  the  subscript  C  refers  to  the  cold  end  of  the  pulse  tube  and 

•  i 

Q  ,     are  the  radiation  and  conduction  losses  associated  with  the 

pulse  tube« 

If  one  now  makes  the  assumption  that  the  period  of  compression, 
T*  ,  is  equal  to  the  period  of  expansion,  T  ,  and,  furthermore,  that 
either  of  these  periods  is  equal  to  vf  T,  where  ^  is  the  percentage  of 
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the  total  time  that  is  devoted  to  compression  or  expansion,  then  equa- 
tions (3,1)  and  (3*2)  can  be  written  as: 


14/ c 


gc    ge. 


W 


+  ^heater  +  \Losses 


-  Qu  Z   0  ,  (3.1a) 


'•H 


and 


y  c 


gc    ge 


+  ^losses  "  QH  =  °  • 


(3.2a) 


By  combining  equations  (3.1a)  and  (3.2a)  one  obtains  the  result: 


QR=  yc_ 


™C  (Tgc-Tge]c-^|T   "T 


gc    ge  w 


(3.3) 


where 


Q  z   Q      + 
R  "  heater 


^ 


losses 


(3.U) 


and   mc    T   -  T   |   is  the  difference  in  enthalpy  flux  at  some 
point  in  the  regenerator  or  pulse  tube.  This  enthalpy  flux  can  be  pre- 
dicted by  Rea's  Theory  as  presented  in  Chapter  II. 

The  difference  in  enthalpy  flux  at  the  cold  end  of  the  pulse  tube 
and  at  the  warm  end  of  the  regenerator  will  now  be  found. 

From  equation  (2.31), 


C  =  2C  [  JL 

1  u-i 


m. 


d£ 

dt 


(3.5) 


and  from  equation  (2.13), 


Cl  =  m  Tgc  "  Tge 


)  • 


(3.6) 
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Now  consider  the  regenerator.  At  the  bottom,  or  cold  end,  x  -   0  and 


•   • 


a  =  hIq  =  n^.  At  the  top,  or  the  warm  end,  x  =  L,  T  =  T^,  and  m  =  in  , 
From  equation  (3«6), 


"V  1  gc    ge 


W 


w 


(3.7) 


R 


Then  from  equation  (3.5>)* 


^^[A 


V 


KAmL 


^C 


R 


dt 


(3.8) 


Note  consider  the  pulse  tube.     At  the  connection  between  the  cooler 

•        •  • 

section  Vtt  and  the  pulse  tube,  x  =  0  and  m  =  m     -  rrijj.     At  the  bottom, 

•       •  • 

or  the  cold  end,  x  =  L,  m  =  m-  =  m^,  and  T  =  T,.  =  Tc#     From  equation  (3.6), 


lTgc-Tge]c  = 


pt 


(3.9) 


Then  from  equation  (3.5)* 


CM       =2C       X     Ljl 
P*  U-l/    \k/ul 


™H 


pt 


dt 


(3.10) 


By  substituting  equations  (3.7)  and  (3.9)  into  equation  (3.3)  one 
obtains : 


0r=  ^cp 


^lptM°i. 


(3.11) 
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Now  by  substituting  equations  (3*8)  and  (3.10)  into  equation  (3.11)  the 
result  is: 


qh  =  2^cp  d: 


dp_ 

dt 


m 


H 


KAyL/  pt  "  "fc  [ka^L 


R 


R 


(3.12) 


Equation  (3.12)  is  identical  with  equation  (3.U)  of  Rea  (h)  with  the 
exception  of  the  ty   factor,  which  apparently  was  omitted  by  Rea.  The 
method  of  determining  this  factor  will  be  commented  on  in  the  following 
chapter. 

Equation  (3*12)  will  now  be  simplified.  One  is  able  to  write  for 
a  perfect  gas: 


pV  =  mRT 


g 


(3.13) 


where  m  is  the  mass.  If  this  equation  is  differentiated  with  respect 
to  time  for  a  particular  volume,  the  result  is  the  general  relation- 
ship: 


™  -  V 
m  -  v 

R 


dp 
dt 


p_   dT^ 
Tg2  dt 


(3.1U) 


Now  consider  the  particular  volume  to  be  V^  the  volume  of  the  cooler 
at  the  end  of  the  pulse  tube.  The  assumption  is  now  made  that  there 
is  perfect  heat  transfer  in  V    This  means  that  the  gas  in  the  volume 
Vtt  is  at  a  constant  temperature  T    With  this  assumption  equation 
(3 .111)  becomes: 
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mH  = 


-  VH 


dp 


(3.1?) 


RT. 


II 


From  equation  (3.13>)  it  follows  that 


m  :  C 
C 


^H  KTK 


dp 
dt 


(3.16) 


By  substituting  equations  (3.15)  and  (3.16)  into  equation  (3*12)  and  by- 
rearranging,  equation  (3.12)  can  be  written: 


QR  =  2^/JL 


dp 
dt 


2Ih/jU     fcpt-ccR' 


MKP!pt 


(3.17) 


where 


r^C     eR       lK£ 


OC   =  Jz 


vt 


"H     €pt      (KPJR 


(3.13) 


and 


8  =   AtLG 


v 


(3.19) 


Equation  (3*17)  is  the  governing  equation  for  the  Pulse  Tubeo 
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CHAPTER  IV 

EXPERIMENTAL  PROCEDURE 
Description  of  Apparatus 

Rea's  Modified  Theory  for  the  Pulse  Tube  was  tested  using  an 
apparatus  constructed  as  shown  schematically  in  Figure  IV.  The  regen- 
erator consisted  of  a  12"-long, stainless-steel  tube,  1"  O.D.  by 
0.020»-wall-thickness,  packed  with  0.0£0»-diameter  lead  shot.  The 
pulse  tube  consisted  of  a  12"-long,  stainless-steel  tube,  5/8"  O.D.  by 
0.016" -wall-thickness.  The  bottom  of  the  regenerator  was  connected  to 
the  bottom  of  the  pulse  tube  by  a  lA"  OJD.  by  0.020" -wall- thickness 
stainless-steel  tube.  Fine-mesh-bronze  screens  were  installed  at  the 
bottom  of  the  regenerator  and  pulse  tube  to  smooth  the  flow.  The  heat 
exchanger  at  the  top  of  the  pulse  tube  was  constructed  from  thin-wall 
l/8"-copper-nickel  tubing.  The  rate  of  flow  of  cooling  water  around 
these  tubes  was  0.5  gal/min.  A  heating  coil  of  high  resistance  wire 
was  wound  around  the  connecting  tubing  between  the  regenerator  and 
pulse  tube  as  a  means  of  introducing  a  known  heat  load  into  the  device. 
Copper-constantan  thermocouples  were  attached  along  the  length  of  the 
regenerator  and  pulse  tube.  Pressure  transducers  were  installed  at 
the  warm  end  of  the  regenerator  and  at  the  water-cooler  end  of  the 
pulse  tube.  The  entire  Pulse  Tube  was  suspended  in  a  vacuum- tight 
chamber  for  thermal  insulation.  During  data  runs  this  vacuum  was 
maintained  at  approximately  2  x  10~5  lfmm   Hg.  Helium  was  used  as  the 
working  gas  for  the  Pulse  Tube.  A  four-speed  motor  allowed  a  choice 
of  cycling  rates  corresponding  to  22.li,  kh*79   82.3,  and  158.0 
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FIGURE  IV 
SCHEMATIC  OF  TEST  APPARATUS 
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reversals  per  minute.  A  complete  description  of  the  apparatus  is 
contained  in  Appendix  D. 

The  important  characteristics  of  the  device  were: 

(li 


(2 
(3 
(li 

(* 

(6 
(7 
(8 
(9 
(10 


VR  =  3.27  in3 
Vpt  =  3.39  in3 
VH    =  0.673  in3 

VptAH  =  $-°h 
W  =  0.96U 
€R  =  0.388 


ept  =  l<>° 


*  =  1.67 
(K^)pt  I  $%   k    ETU/in-hr-ft2-°F 

(K£)R  =  1.1U  x  10*  {^y&        BTU/in-hr-ft2-°F 


(K p)  ^f   assuming  fully  developed  laminar  flow,  was  obtained  from  the 
expression: 


(Kp) 


pt 


-  17.UUK 


where  Nu  =  U.36  and  /3  =    ATL£ 

V 

(K^)q  was  obtained  as  follows:  for  a  spherical  matrix,  P     z     ^    ■  , 

where  d  is  the  diameter  of  the  matrix  spheres.  K  was  determined  from 

R 

Rea's  experimental  results  (U). 
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Method  of  Testing 


Basically  the  objective  of  the  experimental  work  was  to  obtain 
data  from  a  Pulse  Tube  that  could  be  compared  with  the  refrigeration 
effect  which  Rea's  Modified  Theory  predicts  for  that  Pulse  Tube, 

One  of  the  most  important  experimental  results  was  the  axial 
temperature  distribution,  particularly  of  the  pulse  tube.  From  the 
equation  (2.26),  w  is  defined  as: 


w  = 


V 


T^m^R 
o  o 


d£ 

dt 


x 

L 


(lul) 


For  the  pulse  tube,  from  equation  (3«l£) 


Vtt 
m    =    IlVr   S       "  , 


'H 


THR 


dp 
dt 


(U.2) 


Therefore,  for  the  pulse  tube: 


Pt 


0u3) 


vH  h>t 


where  x  is  the  distance  along  the  pulse  tube  measured  from  the  top  as 
shown  in  Figure  III.  For  a  given  pulse  tube,  equation  (U.3)  allows 
one  to  specify  locations  along  its  length  as  a  function  of  w  ^.  This 
was  helpful  in  that  thermocouples  could  be  located  at  any  distance  x 
along  the  pulse  tube  and  that  position  could  be  translated  easily  into 
a  corresponding  w  .  . 

The  assumption  was  made  that  the  temperature  of  the  gas  inside  the 
cooler  was  uniform  and  at  the  same  value  as  that  of  the  cooling  water. 
In  reality  this  is  not  possible,  but  by  using  thin-wall  cooler  tubes 
and  by  allowing  a  short  dwell  time  between  the  compression  and  expansion 
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parts  of  the  cycle,  "perfect"  heat  transfer  could  be  approached.  In 
this  case  (T0)pt  became  simply  TR.  Then  by  measuring  the  temperature 
of  the  pulse  tube  at  each  thermocouple  location  and  by  dividing  each  of 
these  temperatures  by  T^  to  convert  to  the  non-dimensional  temperature 
variable  0  ,  the  experimental  temperature  distribution  of  the  pulse 
tube  could  be  determined  readily.  When  this  experimental  distribution 
was  compared  with  theoretical  temperature  distributions,  one  obtained 
typically: 


6 


C  ef  f 


A 

\\Si- 

i 

C2 

w_ 


w  max 

The  shape  of  all  the  experimental  curves  was  similar  to  the  shape  of 

the  theoretical  curves  except  in  the  vicinity  of  ^mBX,   indicating  that 

the  laminar  flow  assumption  for  the  pulse  tube  was  fairly  accurate. 

In  normal  application  the  theoretical  curves  would  be  entered  knowing 

a  value  of  w^^  and  Qc#  However,  in  this  case,  since  the  entire  pulse 

tube  temperature  distribution  was  known,  it  was  more  appropriate  to  use 

(©c)  ff  as  an  entering  argument.  This  gave  a  more  accurate  value  of 

C   than  if  G„   had  been  usedo 
pt         C 

The  selection  of  the  method  of  comparing  the  calculated  results 

-27- 


with  the  experimental  results  will  now  be  discus sed0  From  equation 
(3oU)  it  should  be  noted  that  QR  includes  the  input  power  from  the 
heater  plus  the  conduction  and  radiation  losses  in  the  regenerator  and 
pulse  tube.  Since  both  conduction  and  radiation  losses  depend  upon  the 
axial  temperature  distribution,  these  losses  would  be  approximately  con- 
stant for  any  operating  conditions  of  the  device  which  resulted  in  the 
same  temperature  distribution.  This  suggested  a  possible  procedure  for 
testing  wherein  a  change  in  the  measured  value  of  QR(Qheater  +<?losses) 
could  be  compared  with  a  change  in  the  calculated  value  of  Q_  for  simi- 

lar  temperature  distributions  and  for  identical  changes  in   ^£   , 

dt 

The  Qiosses>  of  course,  had  to  be  determined.  As  a  first  approxima- 
tion they  were  analytically  obtained  from  equation  (3»17)  for  no 
heater  input.  This  approximation  assumed  that  the  analytical  express- 
ion is  correct o  This  assumption  may  not  be  exact,  but  should  suffice 
as  a  first  approximation.  The  comparisons  suggested  above  could  be 


made  on  a  plot  of  QR  versus 


dt 


2  as  shown  below: 


QTR 


0 


dt 
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A  smooth  curve  cannot  be  legitimately  drawn  through  the  measured  points 


of  the  Qj£  versus 


dp_ 
dt 


2 


plot,  since  there  is  a  third  variable,  C£>  , 
which  was  not  fixed.  Another  method  of  conducting  the  experiment  would 
have  been  to  keep  Q    constant,  and  then  compare  changes  in  experimental 
Qr  with  changes  in  predicted  Q^. 

Based  on  the  above  discussion  an  experimental  procedure  was 
selected  as  follows:  First,  a  cycling  speed  was  arbitrarily  selected 
from  the  four  speeds  available.  Second,  a  ©c  was  arbitrarily  chosen. 
The  selection  of  &     had  one  limitation,  that  being  the  value  selected 
had  to  be  greater  than  the  lowest  possible  ©-,  obtainable  with  the 
device  under  the  limiting  conditions  of  operation  selected.  This  limi- 
tation was  necessary  so  that  a  range  of  operation  was  available  wherein 
the  same  ©c  could  be  maintained  with  an  increasing  heat  load*  Third, 

a  ratio  of  p^/p  .  was  selected  which  would  maintain  the  chosen  ©„. 
max'  *min  C 

Initially  there  was  no  heater  input.  Under  these  conditions  the 
refrigeration  effect  consisted  of  conduction  and  radiation  losses  only. 
It  took  60-90  minutes  for  the  temperature  distribution  to  become  steady. 
Once  steady  state  was  reached,  the 


■2£  required  to  maintain  the 
dt 


selected  C7p  was  noted. 


Next,  the  heater  power  and 


were  increased  in  such  a  manner 


that  the  ©  remained  constant.  It  was  found  that  response  time  to 
slight  changes  in  pressures  was  considerably  less  than  the  response 
time  to  changes  in  the  heater  input.  Therefore,  for  rapid  stabiliza- 
tion, the  heater  power  was  increased  a  fixed  amount  and  then  the  ^£ 

dt 

was  adjusted  to  achieve  the  desired  &n»     Once  heater  power  was  intro- 
duced, it  was  impossible  to  maintain  the  identical  regenerator  tempera- 
ture distribution,  although  ©w  did  remain  nearly  constant.  This,  of 
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course,  infers  that  the  QiOSses  we**e  not  exactly  constant  for  a  series 
of  data  runs.  The  temperature  distribution  of  the  pulse  tube  on  the 
other  hand  could  be  very  closely  duplicated.  Once  the  steady  condition 


was  achieved,  the 


§ 


and  Qneater  were  recorded.  By  adding  the 

0 

Q,  to  the  Qhea+erj  ^e  fy?  was  obtained.     By  continuing  this  process 

a  series  of  experimental  points  of  Qp  versus     dp     c  could  be  obtained, 

rt        <It 

In  the  above  procedure  Qheater  was  determined  by  measuring 
wattage  and  presented  no  difficulty.  Because  of  the  strong . dependency 


of  the  Pulse  Tube  performance  on 


d£ 

dt 


,  provisions  were  made  to  measure 


the  pressure  at  both  the  warm  end  of  the  regenerator  and  water-cooler 
end  of  the  Pulse  Tube.  To  check  the  pressure  drop  through  the  Pulse 
Tube,  the  water-cooler  end  pressure  transducer  was  first  calibrated 
under  static  conditions  using  previously  calibrated  pressure  gauges 
which  were  installed  at  the  warm  end  of  the  regenerator  as  a  reference. 
The  pressure  was  then  allowed  to  cycle.  Photographs  of  the  oscillo- 
scope traces  of  the  pressure  tranducer  output  were  taken  for  different 
pressure  settings  on  the  gauges.  From  these  photographs  and  the 
pressure  transducer  calibration  curve,  the  maximum  and  minimum  pressures 
at  the  water-cooler  end  of  the  pulse  tube  were  obtained.  For  most 
pressure  settings  the  measurements  taken  at  both  ends  of  the  Pulse  Tube 
were  essentially  identical.  For  cycling  speeds  up  to  and  including 
82,3  reversals  per  minute,  the  only  measurable  pressure  drop  for  the 
data  taken  was  approximately  2  psi  at  a  supply  pressure  of  1^0  psigo 
For  these  runs  it  made  no  difference  where  in  the  Pulse  Tube  the 
pressure  measurements  were  taken.  With  the-  highest  speed  of  l£8#0 
reversals,  the  pressure  drop  ranged  from  3  psi  at  £0  psig  to  2k   psi  at 
112  psig.  For  these  runs  the  p^^  and  p  ^^  were  measured  on  the  pulse 
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tube  side  of  the  device o 

Typical  pressure  transducer  traces  are  shown  in  Figure  V.  From 
the  traces  it  can  be  seen  that  the  total  period  can  be  obtained  pre- 
cisely. However,  the  values  for  T       T  ,  and  £[£  are  not  quite  so 

C       e      dt 

obvious*  In  the  analysis  of  the  Pulse  Tube  it  was  assumed  that 

'  c  -  '  e  -  ty^  •  In  keeping  with  this  assumption  the  following  values 

of  yJT   were  decided  upon  by  taking  an  arithmetic  mean  of  T     and  T  : 

c      e 


Reversals/Min 

^f  (sec) 

^ 

22.  k 

1.0 

0.372 

UU.7 

o.$ 

0.372 

82.3 

0.27 

0.372 

158.0 

O.llt 

0.372 

Because  of  the  difference  in  slopes  of  the  compression  and  expan- 
sion parts  of  the  cycle  and  the  slight  non-linearity  of  the  exhaust 


history,  average  values  were  also  used  to  determine  £E 

lar  pressure  history,  Rea  (h)   reported  good  results  by  taking: 


dr 


With  a  simi- 


dp_ 
dt 


Z   ^max  "  ^min 


avg 


£t 


(U.5) 


avg 


where 


*„g  z  vr  =  i  (rc  ,  r 


tt.6) 


and 


p    =  maximum  cycle  pressure 

max  * 


p  .  r  minimum  cycle  pressure, 


A  similar  relationship  for 


d£ 

dt 


was  used  in  this  case. 
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FIGURE  V 


TYPICAL  "pulse  tube"  PRESSURE  -   TIKE  HISTORIES 
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To  determine  a  predicted  value  for  QR,  equation  (3»17)  was  used: 


dp 

dt 


2  V. 


TH     U/? 


£.       /c 


pt 


pt  -   ocCR 


where  oC  is  defined  by  equation  (3.18).  For  any  Pulse  Tube  £     .,    £  R> 
and  V^  would  be  known.  X   and  T  would  also  be  known  depending  on  the 


dt 


were 


gas  used  and  the  temperature  of  the  cooling  water,  ^f  an^ 
determined  by  the  methods  previously  discussed,  (K^?)  ^  was  known 
once  the  pulse  tube  temperature  distribution  became  known.  To  deter- 
mine a  value  of  C  .  a  value  of  the  parameter  <^)  had  to  be  determined. 
From  equation  (2,3£): 


$S2[X]  (-L-]E_i. 

^     2T-1    KA^L  At  T„ 


In  the  above  equation  %   and 


V 


KATL 


were  known;  T0  corresponds  to  TR 


and  was  also  known;  At  was  taken  as  ^T,  which  was  obtained  from 


equation  (U,6);  p  was  taken  as  p .  which  was  defined  as: 

avg 


p„.  =  p  .  +  Pmax  "  Pmin 

*avg   nmn  * 

2 


(U.7) 


Thus  ^>  was  obtained.  Now  since  $,  (w^)  t  and  (^c^eff  were  a11 
known,  the  appropriate  plot  of  0   versus  w  for  the  pulse  tube  was 

entered  and  a  C  ^  extracted.  This  in  turn  led  to  a  value  of  mn/m^ 

•  .0   • 

by  the  procedure  described  in  Chapter  II,  From  mc/mH,  m^,  was  found. 

This  led  to  a  value  for  (K/?)R.  The  next  step  was  to  determine  (wmax)R, 
which  was  obtained  from: 
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(w     )     :   T°  .  Vh  .  -    tt.8) 

mH 


Equation  (1;.8)  follows  from  the  definition  of  wriax,  equation  (2.26), 
assuming  that  m^  is  constant  along  the  connecting  tubing.  All  quantities 

in  equation  (lu8),  then,  were  known.  Now  knowing©,,  and  (w^^)  ,  a 

R 

value  of  CR  was  extracted  from  the  regenerator  curves  of  0  versus  w. 
From  equation  (3»l8)  &*   was  calculated,  since  all  necessary  quantities 
were  known0  Finally,  0  was  calculated  from  equation  (3.17)* 

Sample  Calculation 

Suppose  one  wants  to  predict  the  refrigeration  for  a  Pulse  Tube 
with  the  following  specifications: 

T„  =  cooling  water  temperature  =  £02 °R 
%  =  1|77°R 
V^H  =  0.673  in3 
VR  =  3.27  in3 

Refrigeration  is  desired  at  376°R.  Hence,  TC/TH  z   0»75  and  Tw/Tc  =  1,27, 
Helium  is  to  be  used  as  the  working  gas.  The  cycle  rate  is  to  be  82.3 
reversals  per  minute.  The  pressures  to  be  used  are: 

pmax  =  60  psiS 
Pmin  =  10  Ps±g 
Patm  =  15  psia 
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&   is  calculated  from  equation  (2.35)  to  be  2.1j6.  From  Figures 
XVIII  through  XXVIII  in  Appendix  C  with  w^  =  5.0  and  TC/TH  =  0.75,  a 
cross  curve  of  (^  versus  C  +  can  be  obtained.  Entering  this  cross 
curve  with  (J>  =  2.U6,   a  value  of  C  t  I  0.11*7  is  read.  In  an  analo- 
gous manner  a  cross  curve  of  m  versus  ©  can  be  obtained.  From  this 
curve  mc/mn  =  5.8  is  read.  Then,  from  equation  (U.8), 

(wmax}R  =  1-lfc 

Now  from  Figure  XVI  of  Appendix  C  with  w  =  l.ll;  and  Ty/Tc  Z  1.27, 
CR  =  0.85  is  read.     Suppose  also  that  <X    has  been  calculated  to  be 
0.00811.,  thenoCCR  I  O..OO69.     For  a  cycle  rate  of  82.3,    yT  Z  0.27  sec 
and    y  «  0.372.     From  equation  (U.5): 

J£    =  185  psi/sec  • 

In  addition  (K£)  t  =  Ul.6  BTU/in-hr-ft2-°F. 

Then,  from  equation  (3*17): 

QR  :  h.25  watts 
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CHAPTER  V 

RESULTS 

Presentation  of  Results 

Experimental  data  was  obtained  using  the  apparatus  and  the  pro- 
cedure described  in  the  preceeding  chapter.  Three  series  of  data  runs 
were  made.  In  the  first  series  different  speeds  and  widely  variant 
back  pressures  were  used  for  the  various  runs.  For  all  runs  in  the 
second  and  third  series  a  single  speed  and  a  nearly  constant  back 
pressure  was  used.  Tables  I,  II,  and  III  summarize  the  results  for 
these  series  of  data.  Figures  VI  through  VHI  illustrate  the  compari- 
sons between  the  predicted  and  the  actual  refrigeration  loads.  Observed 
and  theoretical  pulse  tube  axial  temperature  distributions  for  the  three 
series  of  runs  are  presented  in  Figures  DC  and  X.  Figure  XI  shows  a 
plot  of  m  and  C  t  versus  <£  f or  w  s  (■wmax)pt  and  ©  =  (©(Oeff  • 
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Discussion  of  Results 

From  Tables  I,  II,  and  III  it  should  be  noted  that  oCC^  is  less 
than  10/o  of  C  .  in  all  cases  for  this  data  and  for  practical  purposes 
can  be  neglected  when  calculating  QR.  The  data  indicates  that  as  long 
as  the  regenerator  used  in  the  device  is  an  efficient  one,  its  physical 
dimensions  relative  to  the  pulse  tube  dimensions  can  be  smaller  than 
those  of  the  regenerator  used  in  this  experiment  before  its  effect  would 
have  to  be  considered.  If  the  c*CR  term  of  equation  (3.17)  is  neglected, 
that  equation  becomes: 


QR  =  2H>JL 


or,  for  a  given  pulse  tube, 

QR  =  Constant 


op 
dt 


% 


%  \Kfl 


'pt 


(S.D 


pt 


dp 

dt 


pt 


(5.2) 


where  C  ,  is  a  function  of  (j[)  .  All  the  Figures  of  QR  versus 


do 
dt 


in 


this  chapter  are  plotted  using  equation  (5*2).  As  a  result,  all  lines 


of  (p  =  constant  on  the  QR  versus 
ing  through  the  origin. 


dp 


dt 


plots  are  straight  lines  pass- 


For  a  particular  series  of  runs  a  single  temperature  distribution 
was  used  in  determining  a  C   for  any  specific  run  in  that  series, 
since  all  runs  in  each  series  had  practically  identical  pulse  tube 
axial  temperature  distributions.  The  observed  and  theoretical  pulse 
tube  temperature  distributions  are  shown  in  Figures  IX  and  X.  The 
theoretical  results  for  each  series  are  plotted  as  single  curves, 
since  the  deviation  between  runs  in  any  one  series  was  less  than  %%. 
It  is  also  interesting  to  note  that  the  experimental  curves  are  in  good 
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agreement  with  the  theoretical  curves,  indicating  that  the  laminar  flow 
assumption  for  the  pulse  tube  is  a  good  one. 

From  Figures  VI  through  VIII  one  can  determine  approximately  the 
accuracy  of  Rea's  Modified  Theory,  However,  the  accuracy  shown  in  these 
Figures  is  no  better  than  the  accuracy  of  the  variables  involved  in 
arriving  at  these  results.  From  equation  (£.2)  it  can  be  seen  that  QR 
depends  on 


dt 


2 

and  C  .,  The  quantity  C  ,  is,  in  turn,  dependent 


pt°        M"»«w*vj  wpt 


upon  (£   •  In  determining  £?.  and  ^j  ,  average  values  of  pressure  and 
time  were  usedc  The  methods  of  determining  these  average  values  seem 
reasonable.  However,  other  methods  which  would  have  resulted  in  differ- 
ent average  values  could  also  be  justified.  The  importance  of  these 
results,  then,  lies  not  in  attempting  to  show  their  accuracy,  but  in 
showing  that  the  trend  of  the  results  can  be  explained  by  the  theory. 

The  dimensionless  parameters  which  evolve  from  the  theory  are  most 
valuable  in  understanding  the  problem.  For  example,  consider  Runs  2  and 


£  from  Series  1  of  Figure  VI.  Note  that  the 


dp 


dt 


2 

is  approximately  the 


same  for  these  two  runs  and  yet  Qp     for  Run  2  is  more  than  twice 

Al-exper 

Q^     for  Run  5.  This  can  be  explained,  of  course,  by  the  fact  that 
Kexper 

for  Run  $   the  value  of  <j)  was  almost  twice  as  large  as  <|  for  Run  2, 

As  a  result  C  +  for  Run  2  was  more  than  one  and  one-half  times  as  large 

as  C  +  for  Run  £•  An  indication  of  how  C   varies  with  $  in  this  case 

is  shown  in  Figure  XI.  Thus,  Qp     for  Run  5  was  approximately  one- 

"theor 

half  as  large  as  QR     for  Run  2.  The  large  difference  in  <£  for 
these  two  runs  is  attributable  to  the  fact  that  the  back  pressure  for 
Run  5  "was  four  times  larger  than  the  back  pressure  for  Run  2;  also 
Run  5  had  a  \\>T   of  0.27,  while  yT  for  Run.  2  was  0.£0.  The  effect  of 
these  quantities  on  <$  can  be  easily  seen  by  examining  equation  (2.3f>)« 
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Another  example  of  interest  are  Runs  5>  and  k  of  Series  1,  llote 
that  the  total  experimental  refrigeration  for  these  two  runs  is  about 
equal,  yet 


2 

is  more  than  twice  as  large  for  Run  £  as  for  Run  !;• 


dtl 

The  explanation,  of  course,  was  that  (j£   for  Run  5>  was  almost  twice  as 

large  as  <£  for  Run  h.     The  values  of  ^  were  different,  since  <+>T  for 

Run  $  was  0.27  opposed  to  0.5>0  for  Run  i;  and  the  back  pressure  for  Run  $ 

was  20$  greater  for  Run  5  than  for  Run  U.  This  resulted  in  C  .  for  Run  h 

being  more  than  one  and  one-half  times  as  large  as  C  ,  for  Run  5>*  This 

difference  in  C  .  due  to  (§>   was  enough  to  counteract  the  opposite  effect 

on  QPj_,     due  to  <*E 
"theor        dt 

ducing  approximately  the  same  refrigeration  effect. 


p 

and,  hence,  both  conditions  resulted  in  pro- 
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CRA.PT ER  VI 
CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


From  the  analytical  model  of  Chapter  II,  there  appeared  three 
principal  parameters  which  characterize  the  operation  of  a  Pulse  Tube0 
These  parameters  are: 


c  5tf-i  /k^tM   c^  i_ 

^     I  V     1    2m.     dp 

dt 


(2.31) 


w_ 


max 


V 

dp 
dt 

T0moR 

(2.32) 


which  were  described  by  Rea  (h),  plus  a  third  parameter  Q> ,  defined  as 


<£  =  2 


,y-ij  \KA^I  At    t0 


(2.3^) 


The  parameter  <h   resulted  from  the  fact  that  the  ratio 

— 2   for  the  pulse  tube  was  much  closer  to  unity  than 

\ts  it;/  ^  i  y 

for  the  regenerator*  Thus,  the  Simplifying  assumption  which  neglected 
the  term  [p    t^EA    in  the  continuity  equation  (2.3)  in  the  case  of 

the  regenerator  could  not  be  made  in  the  case  of  the  pulse  tube. 

From  the  theoretical  and  the  experimental  results  obtained  from 
the  Pulse  Tube  under  consideration,  it  sras  sppsrent  thst   Pea's  Kc-cifl^-d 
Theory  as  presented  in  Chapter  II  does  describe  the  performance  of  a 
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typical  Pulse  Tube.  There  was  some  deviation  between  the  refrigeration, 
effect  measured  and  the  refrigeration  effect  predicted.  However,  this 
deviation  falls  within  acceptable  limits  for  heat  transfer  data  in  the 
cyclic  case. 

It  was  apparent  from  the  data  obtained  that  there  was  good  agree- 
ment between  the  experimental  and  the  theoretical  pulse  tube  axial  temp- 
erature distributions.  This  leads  one  to  the  conclusion  that  laminar 
flow  in  the  pulse  tube  is  a  good  assumption. 

From  the  analytical  results  obtained  from  the  Pulse  Tube,  it  was 
concluded  that  the  effect  of  the  regenerator  on  the  overall  refrigera- 
tion was  negligible.  Consequently,  the  regenerator  could  have  been 
made  smaller  in  this  case  before  its  effect  would  have  become  signifi- 
cant. 

Experimental  and  theoretical  results  indicated  that  there  is  a 
unique  temperature  distribution  for  a  pulse  tube  operating  with  a 
specified  @q. 

Recommendations 

There  are  several  areas  of  Pulse  Tube  refrigeration  which  should 
be  examined.  First,  the  uniqueness  of  the  axial  temperature  distribu- 
tions of  several  pulse  tubes  with  different  values  of  Wj^^  operating 
with  a  specific  0C  should  be  investigated  over  a  wide  range  of  $  . 
If  this  uniqueness  can  be  demonstrated  for  a  significant  range  of  wmax 
and  $  ,  then  a  suggested  method  for  presenting  the  computer  solutions 
for  the  pulse  tube  differential  equations  is  as  follows: 

(1)  Plot  for  a  fixed  w    a  series  of  cross  curves  of  C_.  versus 
v  max  P& 

$  for  various  values  of  ^q. 

(2)  Introduce  theoretical  axial  temperature  distributions  for 
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some  value  of  <5>  which  covers  the  maximum  feasible  range  in  w. 
VJith  the  temperature  distribution  curves  described  above,  the  0  end 

wmax  °^  any  Pu-^se  'tute  can  ^e  used  to  determine  a  curve  of  constant  C  ^. 
Ey  following  along  this  constant  C  .  curve  to  the  value  of  Wj^^  for  which 
the  cross  curves  of  Cp^  versus  <$)   apply,  a  corresponding  value  of  Op 
can  be  obtained.  This  value  of  Q     represents  a  single  curve  on  the  plot 
of  the  cross  curves  of  Cp^.  and  <|>  .  Thus,  the  relationship  between  Cp^ 
and  <J  becomes  known. 

Second,  a  method  for  determining  an  optimum  size  regenerator  to  be 
used  with  a  given  pulse  tube  should  be  devised. 

Third,  an  optimizing  procedure  for  the  design  of  a  pulse  tube  using 
Rea's  Modified  Theory  should  be  developed. 
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APPENDIX 
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APPENDIX  A 

REA'S  2-PART  MODEL  AND  THE  DEVELOPMENT  OF 
THE  GOVERNING  DIFFERENTIAL  EQUATIONS 

General 

As  a  mathematical  model  for  his  analysis  of  thermal  regenerators 
subjected  to  rapid  pressure  and  flow  cycling,  Rea  devised  his  2-Part 
Model*  From  this  model  Rea  developed  the  differential  equations  that 
apply  to  regenerators  subjected  to  the  above  conditions 0     Rea  proceeded 
to  apply  certain  assumptions  which  simplified  these  equations .  The 
final  result  was  a  system  of  two  differential  equations ,  which  Rea 
solved  numerically.  He  then  reasoned  that  since  a  pulse  tube  is 
simply  an  inefficient  regenerator,  the  same  differential  equations 
should  apply  to  it.  However,  not  all  of  the  assumptions  which  led  to 
the  regenerator  differential  equations  hold  true  in  the  case  of  the 
pulse  tube.  Thus,  the  final  differential  equations  which  govern  the 
pulse  tube  are  not  the  same  as  the  differential  equations  which  govern 
the  regenerator.  This  Appendix  is  broken  into  three  sections.  In 
Section  A.l  Rea's  2-Part  Model  and  his  development  of  the  differential 
equations  for  the  regenerator  will  be  presented  in  detail  as  in  (h)<> 
Section  A  .2  deals  with  the  development  of  the  differential  equations 
for  the  pulse  tube.  The  latter  section  will  not  be  presented  in  as 
rigorous  a  manner  as  Rea's  development  in  Section  A.l,  since  the  two 
developments  are  similar.  All  assumptions  made  in  each  case  will  be 
described  clearly.  Section  A. 3  discusses  the  method  of  solution  for 
the  regenerator  and  pulse  tube  differential  equations. 


A.l 
REA«S  2-PART  MODEL  AND  DEVELOPMENT? 
OF  THE  REGENERATOR  DIFFERENTIAL  EQUATIONS 

Rea's  2-Part  Model 

ReaTs  2-Part  Model,  which  is  shown  in  Figure  II  of  Chapter  I,  is 
defined  as  a  regenerator  in  which  the  mass  flow  rates  at  each  end  of 
the  regenerator  are  in  phase.  By  "in  phase"  it  is  meant  that  during 
each  half-cycle  the  flow  is  always  unidirectional  within  the  regenera- 
tor, i.e.  during  the  pressurization  half-cycle,  mass  flows  in  at  the 
top  of  regenerator  and  exits  at  the  bottom,  and  during  the  expansion 
half-cycle,  mass  flows  in  at  the  bottom  and  exits  at  the  top. 

The  operation  of  Rea's  2-Part  Model  is  as  follows:  the  inlet 
valve  opens  admitting  high  pressure  gas  into  the  regenerator.  The 
pressure  increases  within  the  regenerator  and  the  mass  flow  rate,  m  , 
exits  at  the  bottom.  The  mass  flow  rate  at  the  top,  rn^j,  although 

o 

unknown,  is  certainly  greater  than  mQ  since  some  gas  must  enter  and 
occupy  the  volume  within  the  regenerator.  When  some  maximum  pressure, 

p   ,  is  reached,  the  inlet  valve  closes  and  the  exit  valve  opens. 

max'         ' 

During  the  expansion  half -cycle  mass  flow  equal  to  mQ  enters  at  the 
bottom  and  mass  flow  equal  to  m^  exits  at  the  top.  During  this  half- 

cycle  the  pressure  decreases  to  its  minimum  value,  p  .  •  Once  a 

J  r  mm 

steady,  cyclic  state  is  reached,  the  above  processes  are  repeated  con- 
tinuously with  the  temperature  at  the  top  equal  to  T  and  the  tempera- 
ture at  the  bottom  equal  to  T  . 

o 


It  is.  interesting  to  note  that,  since  m,  is  always  greater  than 

j  top 


m  ,  all  the  gas  entering  at  the  top  does  not  reach  the  bottom.  Thus, 
o' 


it  is  possible  for  some  gas  to  remain  always  within  the  regenerator 
and  flow  back  and  forth  as  the  pressure  cycles. 


FIGURE  XII 


WORKING  SCHEMATIC  OF  THE  REGENERATOR 


max 


Control 
Volume    "- 


Y  mo 


svelopment  of  the  Governing  Differential  Equations  for  the  Regenerator 

The  development  which  follows  with  regard  to  the  regenerator 
differential  equations  is  entirely  attributable  to  Rea  (k)» 

Rea's  initial  assumptions  made  concerning  the  development  of  the 
regenerator  differential  equations  are  as  follows: 
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(1)  The  pressure  drop  in  the  regenerator  is  neglected, 

(2)  The  gas  obeys  the  perfect  gas  lav, 

(3)  The  gas  transport  properties  and  the  matrix  properties  are 
constant, 

(h)     The  cross  sectional  area  of  the  regenerator  open  to  flow,  A  , 
is  constant* 

Figure  XII  shows  a  schematic  of  the  system  to  be  analyzed*  The 

system  consists  of  a  regenerator  of  length  L,  which  is  filled  with  a 

matrix  material  of  large  heat  capacity,  e*g*  spherical  lead  shot*  The 

subscript  m  is  used  to  denote  the  matrix*  The  void  volume  within  the 

regenerator  is  denoted  by  V  and  is  equal  to  A  L*  Curing  each  cycle 

ft  o 

of  operation  the  gas  in  the  regenerator  is  alternately  compressed  and 

expanded  between  the  pressures  p    and  p  .  * 

max     min 

In  accordance  with  the  assumptions  listed  previously,  energy  equa- 
tions will  be  written  for  the  gas  and  for  the  matrix  contained  within 
the  control  volume  of  length  ox,  as  shown  in  Figure  XII, 


Gas  Energy  Equation t 

Energy  flux  into    .  Energy  flux  out  of 
the  Control  Volume    the  Control  Volume 


Rate  of  Change  of 
Energy  within  the 
Control  Volume 


mh 


1-¥'ot§1 


Vog*Vo&($)* 


hjA, 


T  -  T 
g 


dx  + 


f     uA  dx 
'g   © 


or,  upon  simplifying  the  above  equation, 


ilH-Vo^VT^g-*)  ♦10&(V)«0-   <*•« 
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In  equation  (A.l)s  the  first  term  is  the  change  in  enthalpy  flux 
across  the  control  volume,  the  second  term  represents  the  gradient  of 
the  energy  transport  by  gas  conduction,  the  third  term  is  the  convec- 
tive  heat  transfer  between  the  gas  and  the  matrix,  and  the  fourth  term 
is  the  rate  of  change  of  energy  within  the  control  volume. 

Matrix  Energy  Equation ; 


Energy  flux  into 
the  Control  Volume 


Energy  flux  out  of 
the  Control  Volume 


.  Rate  of  change  of 
Energy  within  the 
Control  Volume. 


-k  A  £2  +  h^A 


m  vrs 


+  Mc  $1    dx 
mdt 


or,  upon  simplifying, 


^A!f^A^M  r 


I^AiT  +  ¥T  T  -T   =Mcm  iS 
bx2  \  I  ot 


(A.2) 


In  equation  (A#2)  the  first  term  is  the  gradient  of  the  energy 
transport  by  conduction  along  the  matrix,  the  second  term  is  the  con- 
vective  heat  transfer  from  the  gas  to  the  matrix,  and  the  term  on  the 
right  is  the  rate  of  change  of  energy  of  the  matrix. 

Equations  (A.l)  and  (A.2)  are  energy  balances  for  the  gas  and 
matrix,  respectively.  Subscripts  g  and  m  refer  to  the  gas  and  the 
matrix,  respectively;  h  represents  enthalpy  and  u  represents  internal 
energy;  hm  is  the  heat  transfer  coefficient  (which  has  not  been  assumed 
constant)  and  k  is  thermal  conductivity;  A  is  the  total  cross-sectional 
area  of  the  regenerator  and  AQ   is  the  cross-sectional  area  open  to 
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flow  (AQ  r  €  A) j  A^   is  the  heat  transfer  surface  area  of  the  regenerator 
per  unit  length  (ft2/ft)5  and  M  is  the  matrix  mass  per  unit  length  (lbm/ft). 

The  gas  continuity  equation  will  now  be  developed. 

Gas  Continuity; 

":_-ss  flow  into  the   _   Mass  flow  out  of  the      Rate  of  storage  of 
Control  Volume      "   Control  Volume  mass  within  the 

Control  Volume 


m  =      m+^iUdx      «•  A  \?  Adx 
\  ax        1        it       8° 


• 


am  a  ? 

—  =  -  A     ~~" "" 

ax        on 


lis  .  (A.3) 


Equation  (A.3)  is  the  continuity  equation  and  it  is  perfectly  general 
(within  the  assumptions  listed  previously),  as  it  applies  to  either 
steady  or  unsteady  flow,.  With  the  use  of  the  perfect  gas  law, 


9e    -  -2-    ,  (A.li) 

equation  (A#3)  becomes 

(A.*) 


a  m  _  m  A0   ^ 

ax"    r    3 1 


T 


Equations  (A0l),  (A.2),  and  (A.3)  along  with  appropriate  boundary  con- 
ditions are  the  governing  differential  equations  for  any  type  of  regen- 
erator operation.  Equation  (A.l)  will  now  be  simplified  by  using 
equations  (A0U)  and  (Ac£). 

For  a  perfect  gas,  one  may  write 
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h  =  h0  +  cpTg  (A 06) 

and 


u=h-pv=h-RT  =h  +  c  T     , 
*         g    o    v  g    » 


where  h0  is  the  gas  enthalpy  at  some  reference  temperature • 

The  first  term  in  equation  (A.l)  becomes,  with  the  use  of  equa- 
tions (A.5)  and  (A.6):- 


dxl      I  °ix        P6xl     gl  R        *t\\]  P3x\      g| 

(A.7) 


The  fourth  term  in  equation  (A.l)  becomes,  with  the  use  of  equa- 
tions (A.U)  and  (A.6): 


i  i- 
°«t 


Ll:^  i.  E_\  +  ^d°  >£     .  (A.8) 

S      R   bt\  T„     R   dt 


The  substitution  of  equations  (A.7)  and  (A.8)  into  equation  (A.l) 
yields: 


A0  ilg  +  hfa    T     -  T]     +  fvfo   ^1    =  0       • 
0  by?  \   g         I  R        it 


cpA(mTffj     -kffAQaji  +  h:iAT[Tff.T|     +  ^£0   OP     =  Q       .        (Ao?) 


Equation  (A ,9)  is  the  final  form  of  the  energy  equation  for  the  gasc 
Equations  (A. 9),  (A. 2),  and  (A.5),  along  with  the  appropriate 
boundary  conditions,  completely  describe  the  performance  of  thermal 
regenerators  consistent  with  the  assumptions  listed  on  page  57. 

The  assumption  will  now  be  made  that  terms  involving  thermal 
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conductivity  are  negligible  in  comparison  to  the  other  terms  in  (A09) 
and  (A02)«  That  this  assumption  is  consistent  with  some  practical 
cases  is  shown  by  taking  some  typical  values r 

kg  =  0.07  ETU/hr-ft-°F,   A  =  0o001S2*  ft2,   (t,  -  T)  =  2°F, 

o 

AT  =  U.h3  ftVft,       hT  Z  23^     BTU/hr-ft2-°F,       c    Tg  _  VTo  :  l50°F/ft2. 

b  -J         L2 


Then: 

A2T 
KA0  °     g     "  0.0201*       BTU/hr-ft,   and 

3x2 


hpAy  (T     -  T)  =  2080       BTtJ/hr-fto 

Thus,  the  gas  conduction  term  is  approximately  0.001$  of  the  con- 
vective  heat  transfer  term  of  equation  (A09)o 

In  the  matrix  energy  equation,   (A«2),  with  k^-  20  BTU/hr-ft-°F 
and  A  =  0,005  ft2: 

*2T 

k    A'~   ^    1$     BTU/hr-ft. 

Therefore  the  conduction  term  for  the  matrix  is  less  than  1%  of  the  con- 
vective  heat  transfer  term.  Then  neglecting  axial  conduction,  the  final 
basic  differential  equations  become: 

cp  £  (*Tg)  +  Vt(?8  -  *)  ♦  ^  j|  =  °  (A'10) 


VlK-'   =  *m  H  (A'U) 
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In  equation  (A .10),  the  first  tern  is  the  gradient  of  the  enthalpy- 
flux  along  the  regenerator;  the  second  term  is  the  energy  flux  due  to 
convective  heat  transfer;  the  third  term  is  the  work  of  compression. 

lation  (A „ll)  means  that  the  increase  in  internal  energy  of  the 
matrix  equals  the  energy  given  up  by  the  gas  by  convective  heat  trans- 
fer. Equation  (A ,f>)  is  the  gas  continuity  equation. 

Boundary  conditions  for  the  above  differential  equations  are: 


at  x  Z   0:    m  =  m0(t) 


?g  =  Tg0(t)  (A.12) 

at  x  Z   L:    T  Z   T^j(t) 

Simplification  of  the  Differential  Equations 

Equations  (A. 10),  (Aell),  and  (A.5)  are  three  equations  in  three 

e 

unknowns:  T,  T  ,  and  m.  Thus,  if  p  as  a  function  of  t  is  known,  then 

o 

solutions  consistent  with  the  boundary  conditions,  equations  (A.12), 
will  exist.  Unfortunately,  equations  (A .10),  (A.ll),  and  (A„5)  are  a 
system  of  differential  equations  for  which  no  general  solution  is  known. 
The  method  of  approach,  then,  will  be  to  determine  solutions  which 
satisfy  certain  conditions  over  a  complete  cycle,  but  which  are  not 
exact  in  every  detail  at  each  instant  within  the  cycle 0 

A  cyclic  integral  of  equation  (A.ll),  while  holding  x  constant, 
yields : 
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Q  hAT  j  Tg  -  T]  dt  =  &  Mcm  ^|  dt  =  0.  (A.13) 


5  above  follows  from  the  fact  that  once  steady,  cyclic  conditions 
have  been  established,  the  matrix  temperature  at  a  given  point  in  the 
regenerator  oscillates  around  a  value  of  the  temperature  which  is  a 

tction  of  x  only.  From  this  point  on,  h  will  be  used  for  the  heat 
transfer  coefficient  rather  than  hm. 

A  cyclic  integral  of  equation  (A, 10)  with  x  constant  yields: 


CD 


<fi-(  mTg  j  dt  +    £ hAT  [  Tg  -  TJ     dt  +  ^p     dt  =  0.  (A.lli) 


But,    <T>  <L£     dt  s  0     by  definition,   since  p  is  a  function  of  t  only 

J  a  t 

(no  pressure  drop  in  the  regenerator).  Then  by  using  equation  (A ,13), 
equation  (A. 11;)  becomes: 

Since  mT  will  be  continuous  in  x,  the  order  of  differentiation  and 

to 

integration  can  be  interchanged  and  equation  (A, l£)  becomes: 

i.j(mTg)dt  =  0  , 

or, 

<9    (mTg     dt  =  C0     (for  all  x).  (A .16) 

For  the  2-Part  Model,  the  pressurization  part  of  the  cycle  is 
denoted  by  the  subscript  c  and  the  expansion  part  of  the  cycle  by  the 
subscript  e.  The  mass  flow  rate  at  a  given  point  within  the  regenerator 
will  be  cyclic,  passing  through  zero  at  reversal.  The  mass  flow  versus 
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time  might  appear  as: 


m 


Since  there  can  be  no  net  mass  storage  within  the  regenerator  once 
steady,  cyclic  conditions  prevail, 

(A  ,17) 


f 


mdt  z   0 


and  the  positive  area  must  exactly  cancel  the  negative  area  in  the 
above  diagram.  Similarly,  T^  will  be  cyclic  with  time,  but  in  some 

o 

unspecified  manner  and  <£  T  dt  will  not,  in  general,  equal  zero. 

o 

Then,  the  product  mT  might  have  a  shape  as  follows: 


mT, 


(mTg)e 


where  the  (P  mT  dt  =  C.  from  equation  (A.l6).  This  means  that  the 


i 


positive  and  negative  areas  in  the  above  sketch  do  not,  in  general, 

cancel  each  other.  The  assumption  will  now  be  made  that  mT  versus 

g 

time  can  be  represented  by  the  dashed  lines  in  the  above  figure. 
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Then? 

<CmT  dt  =  (  ^Tgcdtc  ♦   \  neTgedte  =  C£ 


g     'c 


-     ^cTgc^c  +  «eTgere  =  Co     >  (A*l8) 


where  T     is  the  blow  time  during  compression  and  TQ   is  the  blow  time 

during  expansion. 

Furthermore,  if  T  s  T   (compression  period  equals  expansion 

c     e 

period),  equation  (A .18)  can  be  written: 


m  T   +  m  T   =  C.  ,  (A.  19) 

c  gc    e  ge    1         9 


where 

'        cl  =  co/Tc  =  Cc/Te 

From  equation  (A017)> 


f 


mdt  I   \  nigdt  +  \  medt  1  0 


and,  \    mQdt  =     \   medt  .  (A. 20) 

With  the  assumption  that  the  mass  flow  rate  can  be  represented  by  an 
average  flow  rate  which  is  constant  with  time  over  each  half -cycle, 
equation  (A .20)  becomes: 


K  Tc  =  -  k  % 
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and  again  letting  T  z    T      ■> 

c     e 

•    • 

mc  I  me      .  (A. 21) 

Equation  (A. 20)  and,  hence,  equation  (A. 21)  holds  at  any  given  x 
along  the  regenerator.  From  here  on,  from  the  result  of  equation  (A02l), 
m  will  be  used  without  subscript  to  denote  flow  in  either  the  compression 
or  expansion  half-cycle  and  n  is  a  function  of  x  only. 

With  the  assumption  that 

mT  =  (m)(T)   , 
equation  (A'ol°)  along  with  (A02l)  becomes 

■  <Tgc  "  Tge)  =  Cl  <A*22> 

Equation  (A022)  is  a  very  important  and  useful  result.  There  are  no 
restrictions  on  C-,  •  It  can  be  positive,  negative,  or  zero,  depending 
upon  whether  T    is  greater  than,  less  than,  or  equal  to  T_0 
The  usual  heat  transfer  correlation  is  assumed, 

h  =  K.S1  , 

where  n  is  a  positive  exponent  whose  value  depends  on  the  particular 
regenerator  matrix.  In  general  0.f>  -  n  —  0.7.  Hence,  h  is  either 

4 

constant  or  a  function  of  m  only.  Using  equation  (A ,21),  h  will  be  a 
function  of  x  only  and  will  be  written  as  h.  Then  from  equation 
(A.13): 

©  h  A  (T  -  T)  dt  =  h  AT  ®  (T  -  T)  dt  Z  0        , 
since  h  is  independent  of  t,  or, 
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(T  -  T)  dt  =  -  (  (T  -  T)  dt 


(Ao23) 


A  plot  of  (T  -  T)  versus  t  might  appear  as: 

© 

(T  -  T); 


(T  -  T) 


(Tg  -  T)e 


From  equation  (A  .23)*  the  positive  area  on  the  above  sketch  must  equal 
the  negative  area.  As  before,  assuming  (T  -  T)  can  be  represented  by 

o 

some  average  value  and  letting  T    z  T   ,  equation  (A.23)  becomes: 


(T„  -  T)„  =  (T  -  T), 


g 


(A.2ii) 


Since  the  blow-time,  T  9   is  usually  short,  a  reasonable  assump- 
tion to  make  is  that  the  temperature  of  the  matrix  at  a  given  x  will 
vary  very  little  during  a  half-cycle  and,  in  fact,  T  can  be  assumed 
constant  with  time.  In  that  case,  T:  T  (x).  The  validity  of  this 
assumption  can  be  demonstrated  by  the  following  example.  From  equa- 
tion (A.ll): 


h  AT  (^  -  T)  =  Mb-  £T 


6t 


(Aoll) 


For  a  regenerator  12"  long  and  1"  diameter  packed  with  lead  shot  of 
0,050"  diameter,  Mc^L  =  O068  BTU/°F.  Taking  h  =  235  BTU/hr-ft2-°F, 
ATL  =  li.U3  ft2,  and  assuming  (T  -  T)    2°Fj 

o 

hA_  (T    -  T)  z  2080     BOT/hr  =  0.068    $2 

1       &  dt  ' 
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and  jl  =  30,600  Op/hr  =  Q^   op/sec 

th  a  blow  time  of  0.1  second,AT  Z  0o85  °F0  During  a  complete  cycle, 
at  a  given  x,  the  matrix  temperature  would  oscillate  around  c   mean  tem- 
perature with  an  amplitude  on  the  order  of  0.83>  °F,  Thus,  the  assumption 
is  made  that  T  is  independent  of  time  and  is  a  function  of  x  only. 

By  applying  the  above  assumption  to  equation  (Aa2li),  one  is  lead 
to  the  imoortant  result  that: 


T   +  T   =  2T  (A.2£) 

gc    ge 


At  first  glance  equation  (A.25>)  seems  obvious,  but  a  careful  examina- 
tion of  the  preceeding  discussion  will  show  that  two  key  assumptions 
were  necessary  in  deriving  equation  (A025).  The  assumptions  vere 
(l)  the  mass  flow  rate,  m,  is  a  function  of  x  only,  and  (2)  T  is  a 
function  of  x  only. 

By  carrying  out  the  time- differentiation  in  the  continuity 
equation  (A .£),  one  obtains: 


'"-   M1         dJi  .         (A.26) 


£x     R  \  Tg  <*t   Tg fc  at 


The  ratio  of  the  first  term  in  the  brackets  of  equation  (A.26)  to  the 
second  term  is: 


With      AT  -AT  =  0o85  °F   ,  the  above  ratio  becomes: 
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__£.  ~    1.2T   :  MO  to  625. 


ThuSj  the  second  term  in  the  brackets  of  equation  (A026)  is  less  than 
1%   of  the  first  term  and,  consequently,  the  second  term  will  be 
neglected.  Equation  (A„26)  can  now  be  written  as: 


dm 
dx 


RT„  dt 


(A.27) 


Total  derivatives  are  used  in  equations  (A<,27)  since  p  is  a  function  of 
t  only  and  m  is  considered  to  be  a  function  of  x  only.  Furthermore,  T 
will  be  replaced  by  T  in  equation  (A.2?)  since  (T  -  T)  is  small 

o 

(2°  F  to  5°  F).  This  assumption  also  enables  continuity  to  be  satisfied 

over  a  cycle.  Since  m  has  already  been  assumed  to  be  equal  and  opposite 

for  the  compression  and  expansion  half-cycles,  it  is  reasonable  to  expect 

r£  to  be  the  same  for  both  half -cycles.  Then  the  continuity  equation 
cbc 

is  written  as: 


dm 
dx 


RT  dt 


(A„28) 


It  is  now  assumed  that  SE  z     *   constant.  Thus,  pressure  is 

at     ~ 

assumed  to  be  a  sawtooth  function  of  time  as  shown  in  the  sketch  below, 
where  the  zero-axis  is  any  arbitrary  pressure  level: 
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Letting  4?  z     +  constant  also  maintains  consistency  with  the  previous 
assumptions  that  m  and  T  are  functions  of  x  only. 

It  is  now  possible  to  present  the  final  formulation  of  the  regen- 
erator problem.  From  Figure  XII,  the  positive  x  direction  for  the 
regenerator  is  from  bottom  to  top.  Mass  flow  in  this  direction  is  posi- 
tive. During  pressurization  the  flow  is  assumed  to  be  from  top  to 
bottom,  or,  in  the  negati\re  x  direction.  The  gas  energy  equation  \:hich 

was  derived  for  positive  m  with  positive  §E  ,  is  written  for  the  two 

dt  ' 

half-cycles: 


-c 


ps 


mT 


gc 


+  hA, 


Tgc-T 


R 


£p 
dt 


=  0 


(A.29) 


'P  — 

dx 


mT 


i;^ 


+  hAT    Tge-  T 


-  CvA, 


R 


dp 
dt 


Z     0 


(A.30) 


In  equations   (A  .29)  and   (Ao30), 


d£ 
dt 


is  written  as  an  absolute  value 


with  the  plus  or  minus  sign  placed  in  front,  depending  upon  whether 
pressurization  or  expansion  is  occurring  in  the  regenerator.  The  con- 
tinuity relation,  equation  (A .28),  is  written: 


dm 
dx" 


-2 

RT 


dP 
dt 


(A  .31) 


For  convenience,  average  value  bars  have  been  omitted.     Equation   (A. 31) 
holds  for  both  half-cycles,   since  in  the  compression  half-cycle  m  is 
negative,  but    £P     is  positive,  whereas  in  the  expansion  half -cycle 

QX 


m  is  positive,  but     re     is  negative. 


cp 


The  final  formulation  of  the  problem  consists,  then,  of  equations 
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(A«29),  (A.30),  and  (A. 31)  along  with  the  results  previously  attained: 


m  T   -TIC, 

gc    ge      1 


Tgc  +  Tge  =  2T 


(A  .22) 
(A.25) 


Boundary  conditions  for  the  above  set  of  equations  are; 


at  x  :  0:   T  =  T„ 


m  Z   HI 


at  x  =  L:   T  =  TT 


■;] 


(Ao32) 


Equations  (A029)>  (A.30),  (A.31),  (A.22),  and  (A. 25)  are  five 
equations  in  five  unknowns:  m,  T,  T  ,  T  ,  and  C^.  However,  there 
are  only  three  boundary  conditions.  Therefore,  it  is  necessary  to 
eliminate  the  gas  temperatures  from  the  above  equations© 

If  the  differentiation  indicated  in  equations  (A. 29)  and  (A. 30) 
is  carried  out  there  results: 


•cp  i 


•  dTc 


dx 


+  Tgc  — 
dx 


L  OvA-o 
R 

dp 
dt 

=     0 


(A.29a) 


I  m  dTge  +  T  e  £2  J  +  hAT  T   -  T   -  ^o 
dx       dx  /      \       '    R 


dp 

dt 


Z     0   .   (A.30a) 


Now  if  equation  (A#30a)  is  subtracted  from  equation  (A  .2 9a)  the  result 
is: 


-c. 


•   / dT~r       dT„A  dm 

dx  dx  2°         ge/     dx 


+  hAT    Tgc-Tge 


2CyA( 
R 


5 

dt 


=     0        . 


(A.33) 
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Now  by  substituting  equations  (A,22),  (A„25)5  and  (A. 31),  pl^s  the 
relationship  R  z   cp  -  Cy,  equation  (A. 33)  reduces  to: 


*  dT   hAmO,   AA 
m „   T  1  +  __o 

dx   2e  m    c 
P     P 


dp 

dt 


-  0 


(A  .31) 


which  involves  only  the  matrix  temperature,  T. 

The  following  new  variables  are  now  introduced: 


0    -     T/Tr 


m 


-     m/mr 


w     - 


V 


Tjm«R 
0  o 


dp 
dt 


x 
L 


(A.3*) 


Also  the  heat  transfer  coefficient  is  now  assumed  to  be  correlated  by 


h    =     Km11 


(A.36) 


where  K  is  a  constant  with  the  units  BTU/hr-ft  -  jF. 

With  the  aid  of  the  new  variables,  equation  (A.36),  and  the 
definition  ft  z     cp  ,  equation  (A.3I1)  can  be  written: 


2-n  £§L  +  l=k      1-n  - 


m    dw 


m    -  C 


(A.37) 


where 


c  =  y-1  KA^ 


$       2mQV 


dp 
dt 


(A.38) 


The  continuity  relation,  equation  (A„31),  becomes 

dm  -  1 

dw  "  9 


(A.39) 
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Eoundary  conditions  for  equations  (A. 37)  and  (A ,39)  are; 


at  w  :  0: 

m  =  1 

en 

at  w  -  w 

V 

dp 

max 

T0m0R 

dt 

0  =  e. 


w 


(A.liO) 


Equations  (A  .37)  and  (A .39),  along  with  the  boundary  conditions 
UoliO),  are  sufficient  to  determine  m,  T,  and  the  constant  C.  Actually, 
since  everything  in  the  definition  of  C  is  known  except  CL,  determining 
C  is  essentially  determining  C  • 

Thus  the  problem  for  the  regenerator  is  reduced  to  one  of  solving 
a  set  of  two  ordinary  differential  equations  subject  to  the  boundary 
conditions  mentioned.  The  method  of  solution  will  be  discussed  later 
in  section  A. 3. 
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A.2 

THE  GOVFRNT"  |0  :'~-ERr'TIAL  EQUA TIONS  FOR 
REA'S  2-PART  MODEL  AS  APPLIED  TO  A  "pulse  tube" 

A  pulse  tube,  as  long  as  its  heat  transfer  area  per  unit  total 
volume  is  fairly  large,  can  be  considered  to  be  an  inefficient  regen- 
erator. Thus,  a  similar  development  to  that  leading  to  the  final 
differential  equations  for  the  regenerator  could  be  presented  for  the 
pulse  tube.  However,  such  a  development  would  be  repetitious.  The 
procedure  that  will  be  followed  will  be  to  accept  the  assumptions 
listed  on  page  57  and  start  with  equations  (A,°),  (A .5),  and  (A.2), 
which  are  the  gas  energy  equation,  the  gas  continuity  equation,  and 
the  matrix  energy  equation  respectively.  Appropriate  assumptions 
regarding  the  relative  magnitude  of  the  axial  conductivity  will  then 
be  applied.  Next  the  cyclic  conditions  will  be  presented.  At  this 
point  the  relative  magnitude  of  the  terms  of  the  continuity  equation 
will  be  examined.  The  final  formulation  of  the  pulse  tube  problem 
will  then  be  possible 0  The  pulse  tube  differential  equations  will 
evolve  in  much  the  same  form  as  equations  (A. 37)  and  (A .39).  It 
will  be  of  some  assistance  to  present  here  a  schematic  for  the 
2-Part  Model  for  a  pulse  tube.  This  is  shown  in  Figure  XIII.  Note 
that  x  is  measured  positive  from  the  top  down,  which  is  directly 
opposite  from  the  case  of  the  regenerator. 

Equations  (A.9),  (A.£),  and  (A.2)  in  the  form  in  which  they 
were  previously  developed  are  as  follows: 
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c       * 


mT 


"  Vo  ^  *  ^T  (Tg 


-T      +  ^0   ap  =  n 
/  R        ^t 


(A.9) 


<)  m 
3  x 


R      dt 


(A.5) 


FIGURE  Xin 


SCHEMATIC  OF  THE  SYSTEM  TO  BE  ANALYZFD  FOR  THE  "pulse  tube" 
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m. 


Control 
Volume 


max 


fdx 


P  . 
min 


A 


V§*VTJ*.-»--*b{J 


(A. 2) 


As  in  the  case  of  the  regenerator,  the  assumption  will  now  be  made 
that  the  terms  involving  thermal  conductivity  are  negligible  in  compari- 
son with  the  other  terms  in  equations  (A.9)  and  (A.2).  Some  typical 
values  which  substantiate  this  assumption  are  as  f ollows: 

k  =  0.07  BTU/hr-ft-°F,    A^  =  0.0019  ft2,    (T_  -  T)  =  2°F, 
g  o  *     g 
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A     =  0.155  ftVft,         hT  =  6.18     Eru/Wt2-°F, 


|c: —  =  120°F/ft2. 

6  x2  i/ 


Then: 

2 

kgAo  b   Tgr     I     0.016       BTU/hr-ft,       and 


a2 


g 
h^T  (T  -  T)  =  1.92    BTU/hr-ft. 


Thus,  the  gas  conduction  term  is  less  than  1%   of  the  convective  heat 
transfer  terra  of  equation  (A.9)«  In  the  matrix  energy  equation  (A. 2), 
with  k^^:  10  BTU/hr-ft- °F  and  A  =  0.000225  ft2: 

KP   — o  -  °<«270   BTU/hr-ft 
dx2 

the  conduction  term  is  approximately  11$  of  the  convective  term.  If  the 

device  were  operated  with  a  (TH  -  T^  of  190°F,  this  percentage  would 

increase  to  approximately  22#.  Obviously  this  assumption  is  not  as  good 

as  in  the  regenerator  case,  but  without  it  the  problem  would  become 

exceedingly  complex.  Therefore  equations  (A»9)  and  (A.2)  become 

respectively: 

cn  ±    mT  ]  +  hAT  T  -  T]  +  flis  dJL    Z     0  (A.10) 

p  6x  \      s       \  g    /    R   A  t 


and 


The  cyclic  conditions  for  the  pulse  tube  are  identical  with  those 
for  the  regenerator.  Thus,  equation  (A. 22)  applies  equally  as  well  to 
the  pulse  tube  as  to  the  regenerator.  Again  the  heat  transfer 
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correlation  is  assumed  to  be: 

For  the  pulse  tube,  however,  laminar  flow  is  assumed,  fit   the  moment  it 
is  not  intuitively  obvious  that  this  is  a  valid  assumption.  However, 
experimental  data  indicates  that  laminar  flow  is  probably  the  closest 
approximation  that  can  be  made  for  the  pulse  tube.  Thus,  n  =  0,  and 
the  heat  transfer  correlation  for  the  pulse  tube  becomes: 

h  =  constant    • 
Now,  proceeding  with  an  identical  argument  as  that  used  in  the  case 
of  the  regenerator,  one  obtains  equation  (A,2li): 


(Tg  -  T)c  =  -  (Tg  -  T)e  (A.2k) 

Here  the  assumption  is  made  that  the  temperature  of  the  pulse 
tube  wall  at  a  given  value  of  x  will  vary  little  during  a  half -cycle. 
Therefore  T  can  be  assumed  to  be  constant  with  time,  and  T  =  T  (x)  only. 
The  validity  of  this  assumption  can  be  shown  by  the  following  example. 
From  equation  (A .11): 

hAjT  -  T)  =  Me   41  (A .11) 

1  s      m  at 

For  a  stainless  steel  pulse  tube  12"  long  with  O.D.  =  0.625"  and  I.D.  = 
0.591",  Mc^L  =  0.011     BTU/°F,  and  taking  h  =  6.18     BTU/hr-ft2-°F, 
AjL  =  0,155  ft2,  and  assuming  (Tg  -  T)       2°F: 

hA?(T     -  T)  r  1.92       BTU/hr  =  0.011   ||       , 


and 


6  T 

—  =  17U.5  °F/hr.     =     O.OU85  °F/sec. 
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Thus,  with  a  blow-time  of  0,1  second,  A  T  -   0.005°F.  During  a  com- 
plete cycle,  at  a  given  x,  the  pulse  tube  wall  temperature  would 
oscillate  around  a  mean  temperature  with  an  amplitude  of  approximately 
0.005°F. 

Thus,  the  assumption  that  T  is  independent  of  time  and  is  only  a 
function  of  x,  is  a  valid  one. 

By  applying  this  assumption  to  equation  (A.21;),  one  obtains,  once 
again,  the  result: 


Tgc  +  Tge=-2T   . 


(A.2S0 


By  carrying  out  the  time  differentiation  of  the  continuity  equa- 
tion (A .5),  one  obtains: 


d  m 


a0  / 1  d  p  ajg 


(A.26) 


R   T  at   T  c      d  t 
g      g 


The  ratio  of  the  first  term  in  the  brackets  of  equation  (A.26)  to  the 
second  term  is: 

1      6V      'p      <3Tg    ^  ap     Tg 


Tg  ay   Tg     dt 


p    AT 


g 


From  equation  (A922)r 

AT     -  Cl 

g  "  — 

From  equation  (A. 38): 


(A. ia) 


m 


c.  =  2c  J-  ^  -3L 


dp 
dt 


(A«U2) 


By  substituting  equation  (A.U2)  into  equation  (A.U.), 


one  obtains; 
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(A.U3) 


A  Tg  =  2C  _l   !!°    JL     * 

tf-1     m       KATL       dt 

Thus, 

AP     T-     _  1      ^f-1     m       KArpL   At 
p    ATg  "    2C       y      mQ        V         p       8 

Now  using  typical  values  of    tf  =  1.67,     KA-^L/V  =  la.6       BTU/in-hr-ft  -  F, 
C  =  .2,     m/mQ  =  3,  A  t  r  .27  sec,     p  -  75  psia,     Tg  s  390°F,  the  above 
ratio  "becomes: 

A2  h  ~  3 

P  ATg 

For  the  regenerator  this  same  ratio  was  U30  to  62£,  and  the  second  term 
was  neglected.  Rea's  theory  also  neglects  this  term  for  the  pulse  tube. 
However,  the  above  example  indicates  that  the  second  term  cannot  be 
neglected  in  the  case  of  the  pulse  tube.  Equation  (A .26),  then,  remains 
as: 


(A.UO 


Total  derivatives  are  used  because  p  is  assumed  to  be  a  function  of 

t  only,  T  (at  a  given  x)  is  considered  to  be  a  function  of  t  only,  and 

m  is  considered  to  be  a  function  of  x  only.  Also,  since  (T -  T)  is 

D 

small  (2°F0  to  5°F.),  T  will  be  replaced  by  T  in  equation  (A.U3).  Then 
the  continuity  equation  can  be  written  as: 


a™   a0  I  x    dp  _  P  d^ 


dm 

_^o 

L 

dP_ 

P 

dT 

dx 

R 

W 

dt 

V 

dt 

dx    R 


T  dt   T2  dt 
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Now  the  assumption  is  made  that     UJ-f?     can  be  approximated  by 

dt 

.  . S  0  Once  again  4-  is  assumed  to  be  constant  for  both  half-cycles. 
At  dt 

Now  by  substituting  equation  (A.U3)  into  equation  (A.!;£)  and  using  the 

fact  that  4b  -  +  constant,  equation  (A#L£)  becomes: 
at   - 


dm  _  AQ 
dx  "  R~ 


dp 
dt 


-  2C 


2T-1/  4t  T2  U 


or,  simplifying, 


V 
UTLI 


dp 
dt 


(A.U6) 


dm 


-Ao 


dx       RT 


dp 
dt 


1  -  B 


(A.U7) 


where 


B  =  2C  /  j\  P_    1  /  ^o 

^-1/dt     T\ra  /  \  KATL 


(A.U8) 


In  equations  (A.^6)  and  (A.U7)  average  value  bars  have  been  omitted  for 
convenience. 

As  in  the  case  of  the  regenerator,  the  gas  energy  equation  can  be 
written  for  the  two  half -cycles: 


•cp_    ^gc 


+  hAr 


T       -  T 


+  CyA< 


R 


dp 
dt 


=     0 


(A029) 


cp  _     mTge 
dx  \ 


+  hAT     Tge  -  T 


R 


dp 
dt 


Z     0 


(A.30) 


The  final  formulation  of  the  pulse  tube  problem  consists,  then, 
of  equations  (A„U7),   (A. 29),   (A. 30)  along  with  the  results: 
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m  T   -  T   s  C, 
gc    ge    1 


T„_  +  T„a  =  2T 
gc    ge 


(A.22) 
(A.2$) 


Boundary  conditions  for  the  above  set  of  equations  are,  as  for  the 
regenerator: 


at  x  :  0;    T  =  T 


(A.U9) 


m  = 


^o 


at  x  z   L: 


T  =  Tr 


Note  that  the  condition  at  x  :  L  is  different  for  the  pulse  tube  than 
for  the  regenerator  only  in  the  fact  that  the  cold  end  of  the  pulse 
tube  is  at  x  z   L  whereas  the  warm  end  of  the  regenerator  was  at 
x  =  L.  Thus  the  subscript  C  stands  for  "cold"  in  the  pulse  tube  bound- 
ary condition,  whereas  the  subscript  W  stood  for  "warm"  in  the  regen- 
erator boundary  condition,, 

As  in  the  case  for  the  regenerator,  equations  (A.29),  (Ao30), 
(A0U7),  (Ac22),  and  (A.2f>)  are  five  equations  in  five  unknowns: 
m,  T,  T  ,  T  ,  and  C^,  Since  there  are  only  three  boundary  conditions, 
it  is  necessary  to  eliminate  the  gas  temperatures  from  the  above  equa- 
tions • 

If  equation  (A, 30)  is  subtracted  from  equation  (A029)  and  if 
equations  (A.22),  (Ac2£),  and  (A.U7)  are  used  to  eliminate  the  gas 
temperatures,  the  resulting  equation  involves  only  the  temperature  of 
the  pulse  tube  wall: 


;  dT  _  hkfa   +Aq 


dx   2cpm 


dp 
dt 


1-B  |X  -1 

y  -i 


-  o 


(A.SO) 
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Now  if  B  is  replaced  by  equation  (A .I48),  equation  (Ao|?0)  becomes: 


dT 


m 


hA^pC^   Aq 


dx   2cpm    Cp 


dp 
dt 


1  -  2C 


P  1 


k\ 


J-l  At  T 


V 


\  m 


KAtL 


/-l 


*  -1 


-   0  . 

(A.*l) 


The  variables  ©,  m,  and  w,  equation  (A  .35),  plus  the  heat  trans- 
fer correlation,  h  =  K,  are  now  introduced  into  equation  (A„5>l), 
obtaining: 


m2d0+m  U^l  =C  |l  +  *. 


dw 


e 


where 


$  =  2 


/ 


* 


V\  P   1 


(A.52) 


(A .53) 


y-i;  \katl/  At  T 


At  is  Z\t  ,,  for  the  entire  cycle,  and  p  is  p    for  the  cycle  which  is 

avg  '     c  avg        J 

taken  as: 


p    -  p  .  +  pmax  "  min 
avg  "  mm   

2 


The  continuity  equation  (A.U?)  becomes  in  this  case: 


(A.5U) 


dm  .  1  [    <g£ 
dw  "  Q  \       "0m 


(a.55) 


The  boundary  conditions  for  equations  (A0£2)  and  (A05U)  are: 


at  w  Z  0: 

mil 
©=1 

at  w  =  "max 

0  0 

dt 

(A0#) 


«  e=a 
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Thus,  equations  (A .52)  and  (A .55)  together  with  the  boundary 
conditions  (A.56)  describe  the  problem  for  the  pulse  tube. 
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A.3 

METHOD  OF  SOLUTION  FOR  THE  DIFFERENTIAL  EQUATIONS 


Regenerator 


If  equation  (A. 37)  is  solved  for  &&.   , 

dw 

m 


Also, 


¥±    Z    L  •  (A»39) 

dw    © 

The  simplest  method  of  solution  for  the  set  of  equations  (A. 57) 
and  (A .39)  is  to  solve  them  numerically  with  computer  aid.  The 
simplest  approach  is  to  treat  C  as  a  parameter  and  assign  a  range  of 
values  to  it.  The  problem  then  becomes  an  initial  value  problem,  since 
everything  on  the  right  side  of  equations  (A .£7)  and  (A. 39)  at  w  Z   0  is 
known.  The  independent  variable  w  can  then  be  run  out  to  sufficiently 
large  values  to  include  any  desired  w   •  0„  then  becomes  an  unknown 
quantity,  but  can  be  determined  from  a  plot  of Q  versus  w  for  a  given 

value  of  C.  On  the  other  hand,  for  a  given  value  of  ©TT  and  w   ,  the 

'  W     max* 

corresponding  value  of  C  can  be  read  from  the  same  plot. 

Figure  XIV  shows  a  hypothetical  set  of  solutions  corresponding  . 

to  a  family  of  C's.  From  the  figure,  for  C  =  Cp,  O   =  G  at 

w  =  w„„„.  The  entire  regenerator -temperature  distribution  is  deter- 

mined  by  following  the  curve  for  C  z   C_  from  w  =  0  to  w  =  w^«v.  In 

like  manner  one  can  generate  curves  of  m  versus  w  from  equation 
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(A. 39).  In  this  case  the  mass  flow  distribution  can  be  determined  by- 
following  the  appropriate  C  curve  from  w  z   0  to  w  z   w   •  In  addition 

max 

to  specifying  a  range  of  C's,  the  regenerator  matrix  must  be  kn own  as 
well  as  the  working  gas,  so  that  values  of  n  and  2f  can  be  assigned  in 
equation  (A#5>7)» 


FIGURE  XIV 


ILLUSTRATION  OF  COMPUTER  SOLUTTONS  FOR  REGENERATOR 
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max 
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"pulse  tube" 


If  equation  (A. £2)  is  solved  for  —  , 

dw 


C  1  +  -pr  \    -  m 


dS  =  _l_±l_±l_l  (A.S8) 


dw 


m2 


ilso 


dw   ©\    Qm 


Again,  the  simplest  method  of  solution  of  these  equations  is  by 
computer.  Similar  to  the  case  for  the  regenerator,  the  simplest  approach 
is  to  treat  ($)   and  C  as  parameters.  Thus,  for  a  given  (£  =  <^  , 
C  is  assigned  a  range  of  values j  then  a  new  value  is  chosen  for  (J)  and 
C  is  assigned  another  range  of  values.  The  result  is  a  family  of  curves 
of  C  for  a  given  value  of  <§)  ,  as  is  shown  in  Figure  XV.  The  curves  are 
used  in  exactly  the  same  manner  as  that  described  for  the  regenerator. 
To  use  the  curves  one  must  be  able  to  determine  (j>  in  advance.  From 
equation  (A„53)  it  is  apparent  that  to  determine  CJ)  ,  the  following 
information  must  be  known: 

(1)  the  working  gas, 

(2)  the  geometry  of  the  pulse  tube, 

(3)  the  cycling  rate  (See  Chapter  IV ), 
(k)     the  temperature  of  the  cooling  medium, 

(£)  the  average  pressure  at  which  the  device  will  be  operated 
(defined  by  equation  (A.^ii)). 

Details  of  the  computer  program  for  the  regenerator  equations  can 
be  found  in  (ii).  Details  of  the  computer  program  for  the  pulse  tube 
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equations  can  be  found  in  Appendix  B.  Computer  solutions  for  the  regen- 
erator and  pulse  tube  can  be  found  in  .Appendix  C. 


FIGURE  XV 


ILLUSTRATION  OF  COMPUTER  SOLUTIONS  FOR  "pulse  tube" 
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APPENDIX  B 
COMPUTER  PROGRAM  FOR  THE  SOLUTION  OF  THE  "pulse  tube"  EQUATIONS 

General 

Rca  (U)  wrote  a  program  in  Fortran  to  solve  the  regenerator  equa- 
tions (2.29)  and  (2. 30).  A  similar  program  was  written  to  solve  the 
pulse  tube  equations  (2.33)  and  (2.3k).  Runge' s  recurrence  formula  (6) 
was  used  in  the  stepwise  integration 

The  program  was  written  for  ^  s  1.67,  PUMPER)  is  equation  (2.33), 
where  B  corresponds  to  the  mass  flow,  m,  and  R  corresponds  to  the  dimen- 
sionless  temperature  ratio,  ©  •  The  other  symbols  used  in  the  program 
are: 

PAR  =  C 

SPAR  =  $ 

TEMP  =  0  =  T/Tr 

X  -  m  =  m/mc 

SLOPE  =  d  e/dw 

D  =  m  in  Runge  recurrence  formula 

E  s  ©in  Runge  recurrence  formula 

H  =  integration  step  size  =  0.01 

K  =  number  of  integration  steps 

XK  =  number  of  values  of  C  for  each  value  of  $ 

KKK  z     number  of  values  of  ^> 

The  first  data  card  contains  the  initial  values  of  ©  and  m,  the 
step  sise,  the  number  of  integration  steps,  the  number  of  values  of  C 
for  each  value  of  c£>  ,  and  the  number  of  values  of  <$  •  The  remaining 
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•  .* 


data  cards  contain  the  values  of  ^  and  C  for  which  results  are  desired, 


Computer  Program 

*  XEQ 

*  LIBEL 

*  LISTS 

*  FORTRAN 

FUNF  ( B ,  R  )=  ( PAR*  (1 .  +SPAR/R )  -  ( 0  .1<*B ) )/(  B*B  ) 

DIMENSION  A(        ),C(        ),TEMP(        ),Y(        ),SLOPE(        ),W(        ) 

READ  10,TE4P(l),Y(l),H,K,KK,KKK 

DO  25  L=1,KKK 

READ  12, A (L) 

SPAR=A(L) 
3  DO  2£  1=1, KK 

READ  12.0(1) 

PAR=C(I) 
6  DO  20  J=1,K 

DiY(j)+(H*(l./TEMP(j))*(l.-((SPAR*PAR)/(TEMP(J)*Y(j))))/2.) 

E=TEMP(J  )+(H*FUNF(Y(  J),  TEMP(  J )  )/2 . ) 

TEMP(j+l)=TEMP(j)+(H*FUNF(D,E)) 

Y(J+1)=  Y(J)+(R*(1./E)*(1.-((SPAR*PAR)/(E*D)))) 

W(J+l)=H*FLOATF(J) 
20         SLOPE(J+l)=FUNF(Y(J+l),TEMP(J+l)) 

PRINT  1^  H  PAR  SPAR 
2^         PRINT  15'(W(J+1),TEMP(J+1),Y(J+1),SL0PE(J+1),J=20,K,20) 

CALL  EXIT 
10         FORMAT   (3F12.6,IU,I2,I2) 

12  FORMAT  (1F12.6) 

13  FORMAT  (3F12.6)  " 
1$         FORMAT  (IiFl5.5) 

END 

*  DATA 

loO       1.0       .01       600  710 
(OTHER  DATA  FOLLOWS) 
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APPENDIX  C 

THEORETICAL  AXIAL  TEMPERATURE 

AND  AXIAL  MASS  FLOW  DISTRIBUTIONS, 

2-PART  MODEL. 
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FIGURE  XVTT 
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APPENDIX  D 

DESCRIPTION  OF  THE  TEST  APPARATUS 

The  test  apparatus  consisted  of  an  inlet/exhaust  manifold,  Pulse 
Tube,  heater,  heat  exchanger,  and  vacuum  chamber,  A  schematic  of  the 
apparatus  is  shorn  in  Figure  IV  of  Chapter  IV. 

The  inlet/exhaust  manifold  is  shown  in  Figure  XXIX.  The  manifold 
was  constructed  of  standard  -J"  copper-tubing,  using  both  soft  soldered 
and  flared  fitting  connections .  Four  Asco  solenoid-valves  were  used  to 
control  the  flow,  two  on  the  supply  side  and  two  on  the  exhaust  side  of 
the  manifold.  These  valves  had  3/32"  orifice.  The  solenoid  valves  were 
cycled  by  a  cam  operated  micro-switch  arrangement  0  The  cams  were  belt- 
driven  by  a  small,  electric  motor  through  a  four-step  pulley  which 
allowed  cycling  rates  of  22.1*,  Ui.7,  82.3,  and  158.0  reversals  per  minute. 
Two  3/8"  Koke  needle-valves,  one  upstream  of  the  inlet  solenoid-valves 
and  one  downstream  of  the  exit  solenoid-valves,  provided  fine  control 
of  the  gas  flow  rate.  The  check  valve/manual  valve  loops  allowed  for 
reading  maximum  and  minimum  pressures  without  any  pressure  gauge 
fluctuation.  The  maximum  pressure  was  read  on  a  U" -Marsh,  0-300  psi 
gauge.  The  minimum  pressure  was  read  on  a  3"-Marsh,  0-100  psi  gauge. 
The  copper- tubing  of  the  manifold  mated  to  a  1"  0.Do  by  2 J-"-long,  hollow, 
brass  cylinder.  The  lower  end  of  this  cylinder  was  machined  to  allow 
an  "0"  ring  to  be  inserted.  This  end  of  the  cylinder  was  fitted  into 
the  open-end  of  the  regenerator.  Several  short  lengths  of  tubing  were 
fitted  inside  the  brass  cylinder  to  act  as  flow  smoothers. 

The  Pulse  Tube  drawn  to  scale  is  shown  in  Figure  XXX.  It  consisted 
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FIGURE  XXX 
SCALE  DRAWING  OF  PULSE  TUBE 
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of  a  regenerator,  a  void  tube  (pulse  tube),  and  a  short  length  of  connect- 
ing tubing.  The  regenerator  was  constructed  from  a  12n-long,  stainless- 
steel  tube,  1"  O.D.  by  0.020" -vail- thickness  packed  with  G0G£G;,-diameter 
lead  shot.  The  pulse  tube  consisted  of  a  12 "-long,  stainless-steel  tube, 
5/8"  C.D,  by  0.0l6"-vall-thickness.  The  regenerator  and  pul.se  tube  were 
connected  at  the  bottom  by  a  -J"  0,Do  by  0o020" -wall-thickness,  stainless- 
steel  tube.  On  the  pulse  tube  end  of  the  connecting  tubing  was  a  1"  I.D, 
by  l"-long,  stainless-steel  cylinder,  which  x\Tas  filled  with  fine-mesh, 
Maze  screens  to  serve  as  a  flow  smoothing  device 0  Several  layers  of 
this  screening  were  also  placed  at  the  bottom  of  the  regenerator  for  the 
same  purpose,.  The  stainless-steel  components  were  heli-arc  welded  as 
shown  in  Figure  XXX,  The  top  ends  of  the  regenerator  and  pulse  tube 
were  soft  soldered  3"  apart  to  a  brass,  vacuum-jacket,  cover-plate,  which 
was  9"  0oD,  by  3/8"-thick0  Eleven  copper-constantan  thermocouples  were 
installed  along  the  regenerator  and  pulse  tube  at  locations  indicated  in 
Figure  XXX,  All  the  thermocouples  were  led  out  of  the  vacuum  chamber 
through  two  l"-diameter  Latronic  ceramic  seals,  which  were  soldered  into 
pre-driiled  holes  in  the  cover  plate.  The  thermocouple  outputs  were 
recorded  by  a  Leeds  and  Uorthrup  multipoint  recorder©  The  pressure- 
time  history  within  the  Pulse  Tube  was  measured  by  diaphragm- type 
pressure  transducers  mounted  at  the  top  of  the  regenerator  and  pulse 
tube.  An  oscilloscope  was  used  to  indicate  the  output  of  the  pressure 
transducers , 

The  heater  was  constructed  as  follows:  the  %*   stainless-steel 
tubing  connecting  the  regenerator  and  pulse  tube  was  coated  with  a  thin 
layer  of  "Sauereisen  Insa-Lute  Adhesive  Cement",  a  liquid-porcelain  type 
of  material.  After  the  cement  had  hardened,  approximately  two  meters  of 
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a  high  resistance  wire  were  wrapped  on  top  of  ito  The  heater  wire  was 
connected  to  two  32-gauge,  copper,  lead-in  wires „  Then  the  entire 
heater  assembly  was  coated  again  with  a  second  layer  of  cement 0  The 
lead-in  wires  were  lead  through  a  Latronic  ceramic  seal  installed  in 
the  vacuum- jacket  cover  plate.  The  total  resistance  of  the  heater  cir- 
cuit was  121  ohms.  The  heater  circuit  external  to  the  vacuum  jacket 
consisted  of  a  0-lljO  volt  variac  for  fine  voltage  control  0  A  milli- 
ampmeter,  voltmeter,  and  wattmeter  were  installed  in  the  external  cir- 
cuit for  power  measurements. 

The  heat  exchanger  at  the  top  end  of  the  pulse  tube  was  constructed 
as  shown  in  Figure  XXXI.  Thin-wall,  copper-nickel  tubing  1/8"  O.D.  was 
used  to  promote  good  heat  transfer.  The  working  gas  was  funneled  into 
the  tubes,  which  were  approximately  3!""  in  length  and  which  were  dead- 
ended  in  a  brass  tube  sheet  approximately  1"  in  diameter  by  |-f,-thick. 
The  entire  tube  assembly  was  surrounded  by  a  water  jacket,  which  was 
attached  to  the  cover  plate.  Cooling  x^rater  was  pumped  through  the  water 
jacket  at  the  rate  of  0<>5  gallons  per  minute. 
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FIGURE  XXXI 
HEAT  EXCHANGER 
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